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ABSTRACT 

Using our population synthesis code, we found that the typical chirp mass defined 
by (miTO 2 )^/^/(TOi + rfi’2)^l^ of Population III (Pop III) binary black holes (BH-BHs) 
is ^ 30 Mq with the total mass of ^ 60 Mq so that the inspiral chirp signal as 
well as quasi normal mode (QNM) of the merging black hole (BH) are interesting 
targets of KAGRA. The detection rate of the coalescing Pop III BH-BHs is ~I80 
events yr“^(SFRp/(I0“^-^ Mq yr“^ Mpc“^)) • ([fb/(l + fb)]/0.33) -Errsys in our stan¬ 
dard model where SFRp, fb and Errgys are the peak value of the Pop HI star formation 
rate, the binary fraction and the systematic error with Errgys = 1 for our standard 
model, respectively. To evaluate the robustness of chirp mass distribution and the 
range of Errgys, we examine the dependence of the results on the unknown parameters 
and the distribution functions in the population synthesis code. We found that the 
chirp mass has a peak at ^ 30 Mq in most of parameters and distribution functions as 
well as Errsys ranges from 0.046 to 4. Therefore, the detection rate of the coalescing Pop 
HI BH-BHs ranges about 8.3 — 720 events yr“^ (SPRp/(10“^ ® Mq yr“^ Mpc“^)) • 
([fb/(l-l-fb)]/0.33). The minimum rate corresponds to the worst model which we think 
unlikely so that unless (SPRp/(10“^'® Mq yr“^ Mpc“^)) • ([fb/(l + fb)]/0.33) <C O.I, 
we expect the Pop HI BH-BHs merger rate of at least one event per year by KAGRA. 
Nakano, Tanaka & Nakamura (2015) show that if S/N of QNM is larger than 35, we 
can confirm or refute the General Relativity (GR) more than 5 sigma level. In our 
standard model, the detection rate of Pop HI BH-BHs whose S/N is larger than 35 is 
3.2 events yr“^(SPRp/(10“^'® Mq yr“^ Mpc“^)) • ([fb/(l + fb)]/0.33) • Errgys- Thus, 
there is a good chance to check whether GR is correct or not in the strong gravity 
region. 


1 INTRODUCTION 

The second generation gravitational wave detectors such as 
KAGRA\ Advanced LIGO^ Advanced VIRGO® and GEO"* 
are under construction and the first detection of gravita¬ 
tional wave is expected in near future. The most important 
sources of gravitational waves are compact binary mergers 
such as the binary neutron star (NS-NS), the neutron star 
black hole binary (NS-BH), and the BH-BH. As the com- 


^ http://gwcenter.icrr.u-tokyo.ac.jp/en/ 
^ http://www.ligo.caltech.edu/ 

® http://www.ego-gw.it/index.aspx/ 

^ http://www.geo600.org/ 


pact binary radiates gravitational wave and loses the orbital 
energy and the angular momentum, the compact binary co¬ 
alesces. The merger rate of NS-NS can be estimated using 
the binary pulser observation (e.g., Kalogera, Kim, Lorimer 
et al. 2004a,b). However NS-BH and BH-BH merger rates 
cannot be estimated using the observation since no such bi¬ 
naries have been observed so that they can be estimated 
only by the theoretical approach called the population syn¬ 
thesis. For Population I (Pop I) and Population H (Pop II) 
stars, the merger rates of compact binaries are estimated 
by Belczynski, Kalogera & Bulik (2002); Belczynski et al. 
(2007); a (2012); Dominik et al. (2012, 2013). 

In this paper, we focus on Pop HI stars which were 
formed first in the universe with zero metal after the Big 
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Bang. The formation process of Pop III stars has been ar¬ 
gued by many authors such as Omukai & Nishi (1998); 
Bromm, Coppi & Larson (2002); Abel, Bryan & Norman 
(2002); Yoshida, Omukai & Hernquist (2008); Greif et al. 
(2012). The simulations of rotating minihalo in the early 
universe suggest the formation of binaries and multiple star 
systems (Machida et al. 2008; Stacy, Greif & Bromm 2010). 
There are four reasons why we focus on Pop III binaries as 
gravitational wave sources. 

Firstly, Pop III stars are more massive than Pop I stars 
(McKee & Tan 2008; Hosokawa et al. 2011, 2012; Stacy, Greif 
& Bromm 2012) so that Pop III binaries tend to be neutron 
star (NS) and BH. Secondly, the merger timescale of com¬ 
pact binaries due to gravitational wave is in proportion to 
the fourth power of the separation so that even the compact 
binary formed in the early universe can merge at present. 
Therefore, the Pop III compact binary mergers formed in 
the early Universe may be detected by gravitational wave 
detectors such as KAGRA. Thirdly, BH-BH and BH-NS bi¬ 
nary formed from Pop III stars tend to be more massive than 
those formed from Pop I stars while the detectable distance 
of the chirp signal is proportional to 5/6 power of the chirp 
mass defined by (mim 2 )®^®/(mi -|-m 2 )^^® of the compact bi¬ 
nary (Kinugawa et al. 2014). This means that the detectable 
distance of Pop III compact binaries is longer than that of 
Pop I compact binaries so that even if the merger rate per co¬ 
moving volume is smaller, the detection rate can be larger for 
Pop III BH-BH binaries. Fourthly, Pop HI compact binaries 
as gravitational wave source were considered by Belczyn- 
ski, Bulik & Rudak (2004), Kowalska, Bulik & Belczynski 
(2012), they focused on Pop HI stars of mass over hundred 
solar masses or only on the background. However the typ¬ 
ical mass of Pop HI stars is now considered as 10-100 M© 
(Hosokawa et al. 2011, 2012) due to the evaporation of the 
accretion disk by the strong ultra-violet photons from the 
central star. Therefore, Pop HI binary population synthesis 
should be calculated with more realistic lower mass range. 

In our previous study (Kinugawa et al. 2014), we per¬ 
formed Pop HI binary population synthesis with initial mass 
range of 10-100 Mq. The results showed that Pop HI bina¬ 
ries tend to be massive BH ones with chirp mass ~ 30 M©. 
In our standard model, the detection rate of gravitational 
wave by the second generation detectors such as KAGRA is 
140 ■ (SFRp/10-" ® Mq yr-i Mpc"^) • + /6)]/0.33) • 

Errsys events yr“^ where SFRp, fb and Errsys are the peak 
value of the Pop HI star formation rate, the binary frac¬ 
tion® and the possible systematic error due to the assump¬ 
tions in Pop HI Monte Carlo population synthesis, respec¬ 
tively. Errsys = 1 corresponds to our assumption of binary 
parameters and initial distribution functions for our stan¬ 
dard model. However, there are some uncertainties in binary 
evolution and initial distribution function of Pop HI stars so 
that Errsys might not be unity. 


® The definition of the binary fraction is Ni,i„ary/{^single. + 
^binary) where the Nf,i„ary and are the number of sin¬ 

gle stars and the number of binary systems, respectively. Thus, 
SFRp ■ fb/(l + fb) means the binary system formation rate. 


Therefore, in this paper, we perform Pop HI binary pop¬ 
ulation synthesis with several binary parameters and initial 
distribution functions to estimate the variability of merger 
rates and properties of Pop HI compact binaries. That is, 
we estimate the possible range of Errsys to see if the event 
rate is larger than 1 yr“^. In order to perform Pop HI binary 
evolutions, we renewed Hurley’s binary population synthesis 
code (Hurley, Tout & Pols 2002) in our previous paper (Kin¬ 
ugawa et al. 2014). In this paper, we adopt the cosmological 
parameters of (HA,Hm) = (0.6825,0.3175), the Hubble pa¬ 
rameter of Ho = 67.11 km Mpc“^ and the Hubble time 
tHubbie = 13.8 Gyrs (Planck Collaboration 2013). 

The readers of this journal might not be familiar with 
the quasi-normal mode of BH in the title of this paper so 
that we will briefly explain what is the quasi-normal mode 
of BH here. Let us consider the Schwarzschild BH of mass M 
and a small particle of mass m M falling into this BH. If 
the specific angular momentum (L) of the particle is smaller 
than 4GM/c, the particle spirals into the BH. At first the 
gravitational wave depending on L is emitted but in the final 
phase the damped gravitational wave is observed as is shown 
in Figs. 5-4 of (Nakamura, Oohara & Kojima 1987). The 
wave form of this damped gravitational wave is expressed 
by depend only on the mass of 

BH M but not on L. Since ©li > 0 in general, this damped 
oscillation is called the ringing tail or quasi-normal mode. 
Chandrasekahr & Detweiler (1975); Leaver (1985, 1986) ob¬ 
tained the complex frequency of quasi-normal modes includ¬ 
ing Kerr BH case. For the case of Schwarzschild BH, using 
the WKB approximation, Schutz & Will (1985) showed that 
ujr and LOi are determined by the space-time inside ZGM/c?, 
that is, the space-time near the event horizon so that the 
detection of the expected quasi-normal mode of the BH can 
confirm General Relativity in the strong gravity region. If 
it is different from the expected value, the true theory of 
gravity is different from the General Relativity. 

This paper is organized as follows. We describe the mod¬ 
els and the calculation method of Pop HI binary population 
synthesis simulations in §2. Results of our calculation such 
as the properties and the merger rates of Pop HI compact 
binaries are shown in §3. The discussion & summary are 
presented in §4. 


2 THE METHOD OF BINARY POPULATION 
SYNTHESIS SIMULATIONS 

In order to perform Pop HI binary evolutions, we use Pop HI 
binary evolution code (Kinugawa et al. 2014) which is up¬ 
graded from BSE code ® (Hurley, Tout & Pols 2002) for the 
Pop I stars to Pop HI stars case. The result of binary evo¬ 
lutions depends on binary evolution parameters and initial 
distribution functions such as common envelope (GE) pa¬ 
rameters aA, the lose fraction /3 of transferred stellar mass, 
the natal kick velocity, the initial mass function (IMF), the 
initial eccentricity distribution function (lEF) and others 

® http://astronomy.swin.edu.au/ jhurley/ 
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while the details of these parameters and distribution func¬ 
tions will be shown later. We perform Monte Carlo simu¬ 
lations using 10® zero age main sequence (ZAMS) binaries 
for each model. Table 1 shows the parameters and distribu¬ 
tion functions of each model. Each column represents the 
model name, the IMF, the lEF, the natal kick velocity of 
supernova, the CE parameter aA and the lose fraction /3 of 
transfer of stellar matter at the Roche lobe over flow (RLOF) 
in each model. Model name ’’worst” means the worst com¬ 
bination of the parameter and distribution functions each 
model, we calculate 6 cases with different mass range and 
the merger criterion at the CE phase so that the total num¬ 
ber of the models is 84. 

Firstly, we use three mass range cases such as the un- 
derlOO, the overlOO and the 140. In the underlOO case, the 
initial mass range is 10 Mq ^ M ^ 100 Mq. If the stellar 
mass becomes over 100 Mq by the binary interaction, the 
binary evolution calculation is stopped because of no fitting 
formula of the Pop III star evolution for M > 100 Mq (Kin- 
ugawa et al. 2014) due to the lack of numerical data of the 
evolution of Pop III star. This treatment is the same as that 
in Kinugawa et al. (2014). In the overlOO case, the initial 
mass range is also 10 Mq ^ M ^ 100 Mq. However, in 
this case, if the stellar mass becomes larger than IOOMq by 
the binary interaction, the binary evolution calculation con¬ 
tinues by extrapolating the fitting formula of the evolution 
beyond 100 Mq. In the 140 case, the initial mass range is 
10 Mq ^ M ^ 140 Mq. We use the extrapolated fitting 
formula up to M = 140 Mq. The reason for up to 140 is 
that Pop III star of mass larger than 140Mq explodes with 
no remnant (Heger et al. 2003). 

We randomly choose the initial stellar mass from these 
mass ranges with the initial distribution functions such as 
the flat IMF, the log fiat IMF and the Salpeter IMF (See 
§2.1 and 2.2.1). If the stellar mass is over 140 Mq at the 
supernova, the binary evolution is stopped because such a 
high mass star becomes pair instability supernova without 
compact remnant (Heger et al. 2003). Numerical simulations 
of Pop HI star formation by (e.g., Hosokawa et al. 2011, 
2012; Stacy et al. 2012) suggest that the Pop HI protostar 
grows only up to several 10 Mq and the typical mass 
of Pop HI star can be about 40 Mq. In recent simulations 
(Hirano et al. 2014; Susa et al. 2014), the typical mass is 
almost the same as previous study, that is, 40 Mq , however, 
the some Pop HI stars can be more massive than 100 Mq. 
Therefore, we use the initial mass range as 10 Mq ^ M ^ 
140 Mq in the 140 case and study the influences of high 
mass Pop HI binaries for the event rate of gravitational wave 
sources. 

Secondly, we use two merger criteria at the CE phase 
such as the optimistic core-merger criterion and the conser¬ 
vative core-merger criterion (Hurley et al. 2002; Belczynski 
et al. 2002; Kinugawa et al. 2014). In the case of the op¬ 
timistic core-merger criterion, if the condition R'l > 
or R 2 > R'l^ 2 is fulfilled, the primary star merges with sec¬ 
ondary star, where R'l, Rqi, R 2 and Rl _2 a-re the primary 
stellar radius, the Roche lobe radius of the primary star, 
the secondary stellar radius and the Roche lobe radius of 
the secondary star after the CE phase, respectively (Hurley 


et al. 2002). On the other hand, the conservative core-merger 
criterion is R'l + R 2 > at, where at is the separation after 
the CE phase (Belczynski et al. 2002). 

2.1 Brief review of Paper I; our standard model 

In this paper, our standard model is the same as the model 
Hl.f in Kinugawa et al. (2014). In this section, we review 
the model Hl.f in Kinugawa et al. (2014) briefly. The de¬ 
tails are shown in Kinugawa et al. (2014). In order to sim¬ 
ulate the binary evolution, we need to choose initial binary 
parameters such as the primary mass, the mass ratio, the 
separation and the eccentricity. These parameters are de¬ 
cided randomly by the initial distribution function and the 
Monte Carlo method. In our standard model, we adopt the 
flat initial distribution function for the primary mass, the 
flat function for the mass ratio f{q) oc const. (Kobulnicky & 
Fryer 2007; Kobulnicky et al. 2012), the log fiat function for 
the separation f{a) oc 1/a (Abt 1983) and the thermal equi¬ 
librium distribution function of the eccentricity oc e (Heggie 
1975; Duquennoy, Mayor & Halbwachs 1991). We use these 
initial distribution functions except the IMF referenced by 
the Pop I stars observations because there are no obser¬ 
vations suggestions of Pop HI binaries initial distribution 
functions. As for the IMF, some simulations suggest the flat 
or the log flat IMF (Hirano et al. 2014; Susa, Hasegawa & 
Tominaga 2014). Using these initial distribution functions, 
we put the ZAMS binary and start the evolution of the bi¬ 
nary. In order to calculate each stellar evolution, we use the 
formula fitted to the numerical calculations of Pop HI stel¬ 
lar evolutions by Marigo et al. (2001). This fitting formula is 
described by the stellar radius and the core mass as a func¬ 
tion of the stellar mass and the time from the birth of each 
star. The details of fitting equations for Pop HI are shown 
in Kinugawa et al. (2014). We can calculate the evolution 
the binary adding the binary interactions to the Pop HI star 
evolution using the fitting formula. We also need to consider 
the binary interactions such as the tidal friction, the Roche 
lobe over flow, the CE phase, the effect of the supernova 
explosion and the back reaction of the gravitational wave. 

Firstly, we review the tidal friction. The tidal force from 
the companion star changes the stellar radius and the shape. 
In general, the stellar spin angular velocity is different from 
the orbital angular velocity. Therefore the vector of the tidal 
deformation is different from the vector of the tidal force so 
that the tidal torque is generated. The tidal torque trans¬ 
fers the angular momentum from the stellar spin to the or¬ 
bital angular momentum. This interaction changes the bi¬ 
nary separation, eccentricity and the spins of each star. 

Secondly, we review the RLOF. When the star evolves 
and the stellar radius becomes large, the stellar matter is 
captured by the companion star and is transferred to the 
companion star. This phenomenon is called as the RLOF. 
We call the donor star as the primary star. The recipient 
star is called as the secondary star. The region within which 
the stellar material is gravitationally bound to the star is 
called as the Roche lobe radius of that star. If the primary 
stellar radius becomes larger than the Roche lobe radius, 
the primary stellar matter migrates to the secondary star. 
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Table 1. The model parameters. Each column represents the model name, the IMF, the lEF, the natal kick velocity of supernova, the 
CE parameter qA and the lose fraction (3 of transfer of stellar matter at the RLOF in each model. Model name ’’worst” means the worst 
combination of the parameter and distribution functions. 


model 

IMF 

lEF 

natal kick velocity (km s 

aA 

p 

our standard 

flat 

e 

0 

1 

function 

IMFdogflat 

M-i 

e 

0 

1 

function 

IMF:Salpeter 

Salpeter 

e 

0 

1 

function 

lEFxonst. 

flat 

const. 

0 

1 

function 

IEF:e-°® 

flat 

g-0.5 

0 

1 

function 

kick 100 km 

flat 

e 

100 

1 

function 

kick 300 km 

flat 

e 

300 

1 

function 

aX = 0.01 

flat 

e 

0 

0.01 

function 

oA = 0.1 

flat 

e 

0 

0.1 

function 

aX = 10 

flat 

e 

0 

10 

function 

/3 = 0 

flat 

e 

0 

1 

0 

/9 = 0.5 

flat 

e 

0 

1 

0.5 

/3 = 1 

flat 

e 

0 

1 

1 

Worst 

Salpeter 

g-0.5 

300 

0.01 

1 


We define a coefficient /3 as the lose fraction of transferred 
stellar mass descrived as 


M2 = -(l-/3)Mi. 


(1) 


In this case, the accretion rate to the secondary star and 
/3 is determined by the method of Hurley et al. (2002). If 
the secondary star is in the main sequence phase or in the 
He-burning phase, we assume the accretion rate is expressed 
by 

M 2 =-min (10^^, 1) Ml, (2) 

V 'rKH,2 / 


where Mi is the mass loss rate of the primary star and 
is the accretion time scale defined as 


= Ml 

“ Ml ’ 

The Kelvin-Helmholtz timescale rKH .2 is defined as 
GM2(M2 -Mc,2) 


'rKH,2 = 


L 2 R 2 


( 3 ) 


( 4 ) 


where M 2 , Me, 2 , L 2 and R 2 are the secondary stellar mass, 
the core mass, the luminosity and the radius of the star, 
respectively. If the secondary star is in the He-shell burning 
phase, we assume the secondary star can get all matter of 
the primary. Thus, 


iV /2 — —A/i. 


( 5 ) 


In our standard model, we use p function defined by Eq.(2) 
which is computed by the fitting formulae by Hurley et al. 
(2002). However, the accretion rate of the secondary which 
is not a compact object, is not understood well so that we 
use also the accretion rate of the secondary described by the 
constant /3 parameter (For details see §2.2.5). 

If the secondary star is a compact object such as NS and 
BH, we always use /3 = 0 and the maximum of the accretion 
rate is limited by the Eddington mass accretion rate defined 


by 

MBdd = (6) 

= 2.08 X 10-yi + V)-i Meyr-i, 

where kt ~ 0-‘2{l + X) cnl is the Thomson scattering 
opacity and X{— 0.76) is the H-mass fraction for Pop HI 
star. 

When the primary star is a giant and the mass trans¬ 
fer is dynamically unstable, the secondary star plunges into 
the primary envelope and the binary enters the CE phase. 
In this phase, the friction between the primary star and the 
secondary star yields the loss of the angular momentum of 
the secondary to decrease the binary separation, while the 
envelope of the primary is evaporated by the energy liber¬ 
ated through the friction. Consequently, the binary either 
becomes the close binary or merges during the CE phase. 
Now we define ai, af, Mi, Mc,i, Menv,i and Ri as the 
separation before the CE phase, the separation after the 
CE phase, the primary mass, the primary core mass, the 
primary envelope mass and the primary separation, respec¬ 
tively. The separation after the CE phase is calculated by 
the energy formalism (Webbink 1984) defined by 


/ GMc.lMz _ GMiM2 \ _ GMlMenv,! 

\ 2af 2ai J XR-i 


( 7 ) 


where a and A are the efficiency and the binding energy 
parameter, respectively. In our standard model, we adopt 
aA = 1. 

In our calculation, we adopt two merger criteria during 
the CE phase which is shown already in the last paragraph 
of §2 before §2.1. 

When the supernova explosion occurs, the sudden mass 
ejection and the natal kick make the binary obit to change 
drastically. In our standard model, we adopt the natal kick 
velocity equal to zero. In this case, the binary orbit changes 


© 0000 RAS, MNRAS 000, 000-000 



Pop III binary black holes and the quasi normal mode 5 


only by the mass ejection effect in the supernova event. The 
separation and the eccentricity after the supernova explosion 
are described as 


/ 

a 


( |v|^ 

y GMtotal 


GM' 



( 8 ) 



|r X vP 


( 9 ) 


where Mtotai is the total mass, the superscript ' means the 
value after the supernova while v, v and r are the relative 
speed, the relative velocity and the separation vector before 
the supernova, respectively (Blaauw 1961). 

After the binary becomes the compact binary due to 
above binary interactions, the orbit of the binary shrinks by 
the emission of the gravitational waves. The separation and 
the eccentricity are given by 

a GdG^MiMjMtotai 1 + ie" + 
a 5c«a4 (1 - e2)7/2 ’ '■ 


e_ 304G® Ml M 2 Mtotai 1 + 

e ^ 15c5a4 (1 _ e2)5/2 ) 

(Peters & Mathews 1963; Peters 1964). We calculate the bi¬ 
nary evolutions taking account all these binary interactions 
and estimate how many binaries become the compact bina¬ 
ries which merge within the Hubble time. 


2.2 Parameter Surveys 

To estimate the range of Errgys, we calculate the binary 
population synthesis using other initial distribution func¬ 
tions and binary parameters since for Pop III stars we have 
no information on these functions and parameters. However 
in this paper we do not take into account the dependence 
on the initial separation function and the initial mass ra¬ 
tio function because unlike an initial eccentricity functions, 
there are no suggestions for other distribution functions for 
massive binaries although possible dependence of Errgys on 
the change of these initial distribution functions is discussed 
in §4 (Discussion). 


2.2.1 log flat IMF and Salpeter models 

These models correspond to the log flat IMF, that is IMFoc 
dlog(M), and the Salpeter IMF. In recent numerical simu¬ 
lations (Hirano et al. 2014; Susa et al. 2014), IMF of Pop III 
stars might be the log flat IMF. On the other hand, Salpeter 
IMF is acceptable as Pop I IMF. We calculate these IMF 
models in order to estimate the IMF dependence. The other 
initial distribution functions and binary parameters are the 
same as in our standard model. 

2.2.2 lEF'.const. and ® models 

In these models, the IFF are changed from our standard 
model. In general the initial eccentricity distribution might 


be the thermal-equilibrium distribution (IEF:2e) (Heggie 
1975). However in recent observation of massive binaries, the 
eccentricity distribution is not the thermal-equilibrium dis¬ 
tribution. The observation of massive multiple-star systems 
in the Cygnus OB2 (Kobulnicky et al. 2014) implies that the 
observed IFF is consistent with uniform one. On the other 
hand, the observation of massive binaries (M > 15 Mq) 
(Sana et al. 2012) suggests that the power law for the dis¬ 
tribution function of eccentricity as proptoe~^'^. Thus, we 
calculate these two initial eccentricity distribution function 
models. The other initial distribution functions and binary 
parameters are the same as in our standard model. 


2.2.3 kick 100 km s ^ and kick 300 km s ^ models 


The pulsar observations suggest the existence of the NS 
kick. It is observed that the young NSs move with veloc¬ 
ities in range of 200 — 500 km s“^ (e.g. Lyne & Lorimer 
1994; Hansen & Phinney 1997). Since the NS kick veloc¬ 
ity either disrupts the binary or increases the separation, 
the formation rate and the coalescing time of the NS-NS 
and the NS-BH depend on the NS kick velocity. On the 
other hand, the formation rate and the coalescing time of 
the BH-BH might have nothing to do with the natal kick 
velocity because the BH progenitor directly collapses to the 
BH. However, Repetto, Davies & Sigurdsson (2012) suggests 
that the stellar mass BHs have the natal kicks comparable 
to NSs from the distance distribution of the Galactic BH 
(low mass X-ray binaries) above the galactic plane. Pop HI 
BH-BHs are massive, so that they may not have such natal 
kick as stellar mass BHs. However there Is no observation of 
Pop HI BHs. Thus we cannot definitely claim that Pop HI 
BHs do not have natal kicks. Therefore, we take into account 
the natal kick for both NS and BH. In order to estimate the 
dependence on the natal kick, we calculate two models. In 
these models, when stars become compact objects such as 
NS and BH, we assume that the natal kick speed Vk obeys 
a isotropic Maxwellian distribution as 





( 12 ) 


where Uk is the dispersion. In the kick 100 km s“^ model 
and kick 300 km s“^ models, we uses Uk = 100 km s“^ and 
(Jk = 300 km s“^, respectively. The details of the method 
how to calculate the natal kick are shown in Hurley et al. 
(2002). The other initial distribution functions and binary 
parameters are the same as in our standard model. 


2 . 2.4 oA = 0.01, aA = 0.1 and qA = 10 models 

If the primary star becomes a giant and it begins dynami¬ 
cally unstable mass transfer so that the secondary star can 
be engulfed into the envelope of the primary star. In such 
case, the binary enters the CE phase. Once the secondary 
star is swallowed up by the envelope of the primary star, it 
spirals into the core of the primary star due to the orbital 
energy and angular momentum loss by the friction. It is as¬ 
sumed that this spiral-in continues until all the envelope of 
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the primary star is ejected from the binary system. The sep¬ 
aration after the CE phase at is calculated using CE energy 
balance prescription of Eq. 7 

However CE parameters are uncertain. In general, it is 
assumed that a\ — 1 since only a\ is the meaningful param¬ 
eter. Our standard model uses a\ = 1. However, these pa¬ 
rameters can be other values. Therefore, we calculate three 
cases of the CE parameters as aA = 0.01, 0.1, 10. The other 
initial distribution functions and binary parameters are the 
same as in our standard model. 


2.2.5 /3 = 0, /3 = 0.5 and /3 = 1 models 

If the primary star fulfills Roche lobe and it begins dynami¬ 
cally stable mass transfer, the secondary gets the mass from 
the primary star. /3 is called as the lose fraction of trans¬ 
ferred stellar mass defined as Eq. 1. It is considered that 
P varies depending on the binary (e.g. Eggleton 2000). In 
binary population synthesis study, /3 is treated as a func¬ 
tion or the constant parameter. In our standard model, we 
use P as a function (See section 2.1). On the other hand, in 
other studies p is treated as the constant (e.g. Belczynski 
et al. 2002). Thus, we estimate the variabilities of result for 
three constant P cases. The other initial distribution func¬ 
tions and binary parameters are the same as in our standard 
model. Furthermore, if the mass transfer is nonconservative 
{P > 0), the criterion of the stability of the mass trans¬ 
fer should be changed from the criterion of (Kinugawa et 
al. 2014; Hurley et al. 2002) because this criterion assumes 
that the mass transfer is conservative. We use the criterion 
of Eggleton (2011) as 

_ dlogi?L,i 

^ dlogMi 

_ (0.33 -f 0.13gi)(l + qt- pqt) -f (1 - P){ql - 1) - Pqi 


where Rl.i and qt = M 1 /M 2 are the Roche lobe radius and 
the mass ratio. If ((ad = dlogi?ad,i/dlogMi < ((l where Rad 
is the radius when the star reaches hydrostatic equilibrium, 
the binary starts the dynamically unstable mass transfer. 


2.2.6 Worst model 

In this model, we adopt the initial conditions and binary 
parameters which make the worst result in IMF, lEF, kick, 
aA and P (See Section 3). Namely we adopt IMFiSalpeter, 
IEF:e“°'®, kick 300 km s“^, aA = 0.01 and p = 1. We, 
however, think that this worst case is unlikely so that the 
worst case will teach us the minimum merging rate of Pop HI 
BH-BHs. Note that other combination of parameters may 
yield even lower rates. However, we cannot calculate and 
check all3x3x3x4x4 = 432 models. Thus, we choose 
the worst parameters from each parameter region. 


3 RESULTS 


3.1 The properties of Pop III compact binaries 

In order to study the property of Pop HI compact binaries, 
we now show the number of the compact binary formations, 
the number of the compact binaries which merge within 15 
Gyrs and the distribution of the BH-BH chirp mass. The 
details of the difference between each model will be shown 
in the following sub-sections (a) to (g). 

Tables 2 to 15 show the numbers of NS-NS, NS-BH and 
BH-BH binaries for each model from the initial 10® zero 
age main star binary. The meanings of under 100, over 100 
and 140 are explained in the first paragraph of §2. The title 
of each table comes from the change of some parameter or 
that of IMF or lEF from our standard model. In each ta¬ 
ble we also tabulated the number of the compact binaries 
which merge within 15 Gyrs. The numbers in the parenthe¬ 
sis are for the case of the conservative core-merger criterion 
while those without the parenthesis are for the case of the 
optimistic core-merger criterion. The meanings of these two 
criteria were explained in the last paragraph of §2 before 
§ 2.1 

In most cases, the number of merging NS-NSs and merg¬ 
ing NS-BH are very small or zero as one can notice easily 
from Tables 2 to 14. The reasons are as follows. The wind 
mass loss and the mass loss by the binary interactions are 
not so effective for the Pop HI binaries because of the zero 
metallicity and smaller radius of Pop HI stars so that the 
Pop HI binaries tend to disrupt or to increase the separa¬ 
tion by the supernova mass ejection (Kinugawa et al. 2014). 
Therefore we focus on the description of the BH-BHs. Figs. 
1 to 14 show the chirp mass distribution of BH-BH binaries 
which merge within 15Gyr for each model. In each hgure, the 
red, green, blue, pink, light blue and grey lines correspond to 
underlOO case with optimistic core-merger criterion, overlOO 
case with optimistic core-merger criterion, 140 case with op¬ 
timistic core-merger criterion, underlOO case with conserva¬ 
tive core-merger criterion, overlOO case with conservative 
core-merger criterion and 140 case with conservative core¬ 
merger criterion, respectively. One can see that in all models, 
the peak of the observable chirp mass distribution is about 
30 Mq. Pop HI binaries with each mass M < 50 Mq are 
unlikely to be the GE phase. They evolve via some mass 
transfer phases and their mass loss is smaller than the evo¬ 
lution passes via a GE phase. They tend to lose 1/10-1/3 of 
their mass so that they tend to be 20-30 Mq BH-BHs. Pop 
HI binaries with each mass M > 50 Mq are likely to be the 
GE phase and they lose 1/2-2/3 of their mass so that they 
tend to be 25-30 Mq BH-BHs too. Therefore, the peak of 
chirp mass become 25-30 Mq. 

Figs. 15 to 20 show the dependence of merger time dis¬ 
tribution of BH-BH binaries for each model. In each figure, 
we describe underlOO cases with optimistic core-merger cri¬ 
terion because the merger time distributions do not change 
a lot in other cases and core-merger criteria. The most im¬ 
portant characteristic is that the merger rate for t > 10®'®yr 
is almost constant for every model. 
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(a)our standard model 

UnderlOO case is the same as the result of the model Ill.f 
in Kinugawa et al. (2014). OverlOO case is equal to the re¬ 
sult of underlOO case plus the binaries whose star becomes 
more massive than 100 Mq by the mass transfer so that 
the number of massive BH-BH mergers in overlOO case is 
more than that of the underlOO case (See Fig. 1). Thus, the 
number of Pop III BH-BHs which merge within 15 Gyrs in 
overlOO case increases about ten percent compared to the 
that of underlOO case (See Table 2). However, the peak of 
the chirp mass is not changed (See Fig. 1) In 140 case the 
peak of the chirp mass distribution (Fig. 1) is almost the 
same as that of underlOO and overlOO cases, but the chirp 
mass distribution of 140 case has tail in the region of 25-60 
Mq since the high mass Pop III binaries (>100 Mq) tend 
to be high mass BH-BHs (40-60 Mq) via the CE phase. The 
event rate of high chirp mass BH-BHs is large because the 
detectable volume (V) is H oc The tail of chirp mass 

distribution dN/dM in 140 case is proportional to in 

the range of 25-60 Mq. 


(b) log flat and Salpeter IMF models 

In these models, the initial masses tend to be smaller than 
that of our standard model. Thus, the numbers of the BH- 
BH formations and merging BH-BHs decrease because the 
BH-BH progenitors decrease due to IMFs (See Table 3 and 
4). On the properties of the chirp mass distributions in the 
log flat and Salpeter IMF models, the number of merging 
BH-BHs where each mass is more massive than 30 Mq is 
smaller than our standard model due to the steepness of 
IMFs. However, the peak mass is independent on the IMF 
(See Figs. 2 and 3). In 140 cases, the chirp mass distribution 
dN/dM of IMF for log flat model is proportional to M~^ 
and that of Salpeter model is proportional to ® in a 
range of 25-60 Mq. In Fig. 15, the shape of each merger 
time distribution is almost the same independent of IMFs, 
although the number of BH-BHs decreases dependent on the 
steepness of IMF. 


(c)IEF:const. and e models 

In these models, the initial eccentricities tend to be smaller 
than that of our standard model. If the initial eccentricity is 
small, the decrease of the separation by binary interactions 
such as the tidal friction and the gravitational radiation is 
suppressed. Thus, the merger rate by the binary interaction 
decreases. Therefore, the number of BH-BH formation rate 
increases while the number of merging BH-BHs deceases al¬ 
though the influence to the merger rate is very small (See 
Table 5 and 6). The properties of the chirp mass distribu¬ 
tions are almost the same as in our standard model (See Fig. 
4, 5). The properties of the merger time distributions are al¬ 
most the same as in our standard model too (See Fig. 16). 
We can say that the effect of the eccentricity distribution is 
not so large. 


(d)kiek 100 km s ^ and kiek 300 km s ^ models 

In these models, the natal kick disrupts binaries or makes 
their orbits wide and eccentric. Firstly, we argue the kick 
100 km s“^ model. The number of the BH-BH formations 
decreases by the natal kick compared to our standard model. 
However, the number of merging BH-BHs becomes about 1.5 
times larger than that of our standard model (See Table 7). 
The merger timescale by the gravitational radiation is given 
by. 

^ 5 /GMtotai 

W(eo = 0) = — - j j 

= inio f VV yY ^totai V^ 

^ V43R0y^ VSOMq^ VsOMQy) VoOMQy) 

(14) 


where ao and eo are the separation and the eccentricity when 
the BH-BH is formed. Thus, the BH-BHs whose separation 
is larger than about 50 R© cannot merge within Hubble 
time. On the other hand, the escape velocity and the orbital 
velocity of the binary system are given by 


"t^esc 


2GM 


= 500 km s"^ 


/ M 
V30 MqJ 


a 

43 R© 


- 1/2 


(15) 

(16) 





= 350 km s ^ 


M 

30 Mq 


1/2 


43 Rf 


0 


- 1/2 


(17) 

(18) 


where a is the separation. Equations 15-18 tells us that even 
though the kick velocity which is 100 km s“^ is aligned 
with the orbital velocity, the BH-BH progenitors whose sep¬ 
aration is smaller than 43 R© cannot be disrupted by the 
natal kick. Therefore, the BH-BHs disrupted by the natal 
kick of 100 km s~^ rarely contribute to the merger rate of 
Pop HI BH-BHs from the beginning. While the increase 
of the eccentricity by the natal kick reduces the merger 
timescale due to the gravitational waves since Emerge (e) ~ 
(1 — e^)^^^rmerge(eo = 0) (Peters & Mathews 1963; Peters 
1964), so the number of merging BH-BHs is larger than 
that of our standard model. The chirp mass distribution 
is, however, almost the same as that of our standard model 
(See Fig. 6). In Fig. 17, the merger time distribution of the 
kick 100 km s“^ model is almost the same as in our stan¬ 
dard model. However, the number of merging BH-BHs which 
merges at the early universe is larger than that of our stan¬ 
dard model due to the increase of the eccentricity by the 
natal kick. 

Secondly, we argue the kick 300 km s“^ model. The 
number of the BH-BH formations decreases by the natal kick 
compared to our standard model and the kick 100 km s“^ 
model. The number of the merging BH-BHs also decreases 
unlike the kick 100 km s“^ model. The reasons are that 
the sum of the orbital velocity and the kick velocity some¬ 
times can exceed 500 km s“^ and that if the binary was not 
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disrupted, the separation of compact binaries tends to ex¬ 
tremely widened by the high natal kick velocity. Thus, the 
many binaries will be disrupted and the coalescence time 
of the survived binaries due to the emission of the gravi¬ 
tational wave tends to be so long that the Pop III binaries 
cannot merge within 15 Gyrs. Especially, the low mass bina¬ 
ries are more susceptible to this effect than the massive bi¬ 
naries. The chirp mass distributions is also changed and the 
peak mass region moves to a little more massive than that 
of our standard model (See Fig. 7). In Fig. 17, the merger 
time distribution of the kick 300 km s“^ model is different 
compared to our standard model. The number of BH-BHs 
which merge at the early universe is larger than that of our 
standard model due to the increase of the eccentricity by 
the natal kick. On the other hand, the progenitors of BH- 
BHs which have large separation but can merge within the 
Hubble time in our standard model tend to be disrupted by 
the natal kick. Thus, the number of merging BH-BHs which 
have large merger time decrease. 

(e)aX = 0.01, aX = 0.1 and aX = 10 models 

In these models, the CE phase results in the merger of binary 
or the change of the separation. The separation after the 
CE phase is determined by the CE parameter aX. If aX is 
small, the separation after the CE phase tends to be small 
and to merge during the CE phase due to the increase in the 
loss of the orbital energy. Easily the binary does not survive 
during the CE phase. But if binary survives during the CE 
phase, the binary easily merges by the large luminosity of 
GW emission due to the tight orbit. On the other hand, if aX 
is large, the separation after the CE phase tends to be large 
due to the decrease in the loss of the orbital energy. Thus, 
the binary tends to survive during the CE phase. When the 
binary survives in the CE phase, however, it is hard to merge 
due to the small luminosity of GW. 

Firstly, we consider aX = 0.01 model. In this case, the 
number of the BH-BH formations decreases and the num¬ 
ber of the merging BH-BHs becomes about 1/3 as large as 
our standard model. The parameter aA is so small that the 
almost all binaries merge during the CE phase. The Pop 
HI giant with initial mass larger than 50 M© can enter the 
CE phase and these binaries tend to merge. Thus, the sur¬ 
vived binaries evolved via the RLOF. The binaries which 
evolved via the RLOF tend to be smaller than 50 Mq and 
the massive binaries merge more easily. Thus, in the chirp 
mass distribution (Fig. 8), the number of the merging high 
mass BH-BHs becomes smaller than our standard model. 
Especially, 140 case is easily affected by these effects. In 
Fig. 18, the shortest merger time is larger than that of our 
standard model, because almost all progenitors of BH-BHs 
merge during the CE phase. 

Secondly, we consider the aA = 0.1 model. In this case, 
the number of the BH-BH formation decreases compared 
with our standard model. While, the number of the merging 
BH-BHs is almost the same as that of our standard model. 
Owing to the small aA, the binary does not tend to sur¬ 
vive during the CE phase. But if binary survives in the CE 
phase, the binary becomes close binary and easy to merge 


due to the emission of GW. Thus, the number of the merg¬ 
ing BH-BHs does not decrease. In 140 case, the number of 
merging BH-BHs which are low mass becomes larger than 
our standard model (See Fig. 9) because the low mass BH- 
BHs can merge more easily than our standard model due to 
moderately small aA. In Fig. 18, the shortest merger time is 
smaller than that of our standard model due to moderately 
small aA. aA is small but the binaries which become the CE 
phase do not tend to merge during the CE phase and they 
have close orbit due to small aA. 

Thirdly, we consider the aA = 10 model. In this case, 
the number of the BH-BH formation increases compared 
with our standard model. While the number of the merging 
BH-BHs decreases compared with our standard model. aA 
is so large that the separation after the CE phase tends to 
become large and the binary does not merge during the CE 
phase. Thus, the number of the BH-BH formations increases. 
But, due to the large separation the BH-BH binary does not 
merge within 15 Gyrs. Therefore, the number of the merg¬ 
ing BH-BHs decreases. On the other hand, the number of 
merging high chirp mass BH-BHs is more than our standard 
model (See Fig. 10), because the BH-BHs which are formed 
after the CE phase tend to have wide orbit due to high aA 
and the wide massive BH-BHs can merge more easily than 
the wide low mass BH-BHs. Especially, in 140 case this ef¬ 
fect is remarkably clear. In Fig. 18, the shortest merger 
time becomes larger than that of our standard model and 
the number of merging BH-BHs which merge at the early 
universe is smaller than that of our standard model due to 
the same reason. 

(f)/3 = 0, /3 = 0.5 and /3 = 1 models 

In these cases, the accretion rate during the RLOF is 
changed by the loss fraction of transferred stellar mass f}. 
Firstly, we consider the /3 = 0 model. 

The result of our standard model and the result of /3 = 0 
model are the same (Table 12,2 and Figs. 1, 11, 19). We use 
the Hurley’s fitting p which is fitted by Pop I binaries for 
the Pop HI case. We found the Hurley’s fitting /3 is same as 
/3 = 0 in Pop HI case. Fig. 11 and Fig.l are the same. This 
means that the RLOF of our standard model (equation 2) 
is the conservative mass transfer prescription. 

Secondly, we consider /3 = 0.5 model. In this case, the 
mass loss during the mass transfer makes the separation 
wide so that the mass transfer tends to be dynamically sta¬ 
ble. The number of the binaries which are merged during 
the CE phase decreases. Therefore, the number of BH-BH 
formation increases compared to our standard model. How¬ 
ever, the number of merging BH-BHs decreases because the 
binary does not tend to be a close binary by the CE phase. 
The peaks of the chirp mass distributions are almost the 
same as in our standard model (See Fig. 12). While the high¬ 
est mass of merging BH-BHs decreases compared with our 
standard model because of the mass loss during the RLOF. 

Thirdly, we consider /3 = 1 model. In this model, the 
number of BH-BH formation is larger than our standard 
model, but it is smaller than that of /3 = 0.5 model. Like 
P = 0.5 model, the mass transfer tends to be dynamically 
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stable. Especially, in this case, there are no paths to the 
CE phase via dynamically unstable mass transfer. There are 
only paths to the CE phase in which the secondary plunges 
into the primary due to the eccentric orbit or in which the 
each star is a giant and the binary become the contact bi¬ 
nary due to the expanding. Thus, the number of the bina¬ 
ries which are merged during the CE phase decreases, so 
that the number of BH-BH formations increases compared 
with our standard model. However, the mass loss during the 
mass transfer is so large that some stars which can be BH 
by the mass accretion during the RLOF cannot be a BH 
but a NS. Therefore, the the number of BH-BH formations 
decreases compared with /3 = 0.5 model, but the number of 
NS-BH formations increases accordingly. Furthermore, since 
the progenitors of merging BH-BHs are hard to enter the CE 
phase and tend to become wide orbits due to the mass loss 
during the mass transfer. Thus, the number of merging BH- 
BHs decrease. The major progenitors of merging BH-BHs do 
not enter the CE phase and lose their mass during the mass 
transfer. But, since the Pop HI statr radius is small (See 
Fig.2 in Marigo et al. (2001) and Fig.l in Kinugawa et al. 

(2014)), the mass loss during the RLOF tends to stop right 
away and the separation tends to be close enough that Pop 
HI BH-BHs can merge within 15 Gyrs. The highest mass 
peak region of the chirp mass distributions becomes smaller 
than our standard model and the highest mass of merging 
BH-BHs decreases due to the mass loss during the RLOF 
(See Fig. 13). 


(g) Worst model 

In this model, we choose the initial conditions and binary pa¬ 
rameters which will make the worst result in (b)IMF, (c)IEF, 
(d)kick, (e)aA and (f)/9. Thus, we adopt (b)IMF:Salpeter, 
(c)IEF:e“° ®, (d)kick 300 km s“^, (e)aA = 0.01 and (f)/3 = 
1. Especially, we already know that (b)IMF, (d)kick, (e)aA 
and (f)/3 influence the result very much so that the result 
of BH-BH formation and the number of merging BH-BHs 
are determined by these parameters and IMF. Each effect 
of (b)IMF:Salpeter, (d)kick 300 km s“^, (e)aA = 0.01 and 
(f)/3 = 1 makes the number of merging BH-BHs decrease 
(See Table 4, 8 and 14). Thus, the number of merging BH- 
BHs extremely decreases compared with our standard model 
(See Tabel 2 and 15). On the properties of chirp mass distri¬ 
bution, the number of merging BH-BHs more massive than 
30 Mq decreases and the gradient of the chirp mass dis¬ 
tribution of merging BH-BHs is much steeper than that 
of our standard model (See Fig. 14 (a) (b)). The Salpeter 
IMF makes the number of high mass stars to decrease and 
non-conservative mass transfer (/? = 1) prevents to merge 
BH-BHs while in our standard mode they can be merging 
BH-BHs after the CE phase. Furthermore, even though the 
massive binaries become the CE phase, they usually merge 
during the CE phase due to the very small a\. Therefore, 
the peak region of the chirp mass is 25-30 Mq even in 140 
case. 
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Table 2. our standard model 

This table shows the numbers of NS-NS, NS-BH and BH-BH binaries and the numbers of each compact binary which merges within 15 
Gyrs for our standard model. 15Gyrs is used in order to compare our results with previous works. To estimate the present merger rates, 
we use 13.8Gyrs as the present age of the universe. The meanings of under 100, over 100 and 140 are explained in the first paragraph of 
§2. The numbers in the parenthesis are for the case of the conservative core-merger criterion while those without the parenthesis are for 
the case of the optimistic core-merger criterion. The meanings of these two criteria were explained in the last paragraph of §2 before §2.1 


Table 3. IMFdogfiat 


Table 4. IMF:Salpeter 


Table 5. IFF:const. 



underlOO 

overlOO 

140 

NS-NS 

0 (279) 

0 (279) 

0 (195) 

NS-BH 

185335 (187638) 

185335 (187638) 

153435 (155694) 

BH-BH 

517067 (522581) 

534693 (540316) 

595894 (604930) 

merging NS-NS 

0 (279) 

0 (279) 

0 (195) 

merging NS-BH 

50 (149) 

50 (149) 

825 (1255) 

merging BH-BH 

115056 (120532) 

131060 (136645) 

128894 (137903) 


Same as Table 2 but for IMFilogflat model. 


underlOO 

overlOO 

140 

NS-NS 

NS-BH 

BH-BH 

merging NS-NS 
merging NS-BH 
merging BH-BH 

2 (789) 

168100 (169794) 
350169 (353524) 

2 (789) 

68 (183) 

74745 (78054) 

2 (789) 

168100 (169794) 
357989 (361378) 

2 (789) 

68 (183) 

81786 (85129) 

1 (693) 

157106 (158831) 
405922 (410802) 

1 (693) 

374 (579) 

87590 (92450) 


Same as Table 2 but for IMF:Salpeter model. 


underlOO 

overlOO 

140 

NS-NS 

NS-BH 

BH-BH 

merging NS-NS 
merging NS-BH 
merging BH-BH 

5 (1994) 

93085 (93793) 
132534 (133485) 

5 (1994) 

64 (164) 

25536 (26468) 

5 (1994) 

93085 (93793) 
133880 (134835) 

5 (1994) 

64 (164) 

26720 (27656) 

3 (1957) 

92861 (93603) 
144096 (145294) 

3 (1957) 

97 (216) 

28378 (29564) 


Same as Table 2 but for IMFdogfiat model. 


underlOO 

overlOO 

140 

NS-NS 

0 (358) 

0 (358) 

0 (255) 

NS-BH 

183460 (184761) 

183460 (184761) 

152099 (153548) 

BH-BH 

522809 (526892) 

541264 (545459) 

602071 (608210) 

merging NS-NS 

0 (358) 

0 (358) 

0 (255) 

merging NS-BH 

43 (130) 

43 (130) 

843 (1087) 

merging BH-BH 

111106 (1115171) 

127904 (132081) 

124714 (130831) 
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Table 6. IEF:e 


Same as Table 2 but for IEF:e model. 



underlOO 

overlOO 

140 

NS-NS 

0 (365) 

0 (365) 

0 (258) 

NS-BH 

181650 (182388) 

181650 (182388) 

150779 (151805) 

BH-BH 

523285 (526534) 

542015 (545389) 

602575 (607054) 

merging NS-NS 

0 (365) 

0 (365) 

0 (258) 

merging NS-BH 

38 (100) 

38 (100) 

774 (964) 

merging BH-BH 

107594 (110832) 

124620 (127983) 

121494 (125955) 


Table 7. kick 100 km s ^ 


Same as Table 2 but for kick 100 km s ^ model. 



underlOO 

overlOO 

140 

NS-NS 

283 (794) 

283 (794) 

180 (516) 

NS-BH 

32701 (34778) 

32701 (34778) 

32014 (34144) 

BH-BH 

191755 (197327) 

208268 (213962) 

234117 (243348) 

merging NS-NS 

17 (526) 

17 (526) 

6 (342) 

merging NS-BH 

2527 (3016) 

2527 (3016) 

3218 (3762) 

merging BH-BH 

117415 (122830) 

132066 (137603) 

135758 (144554) 


Table 8. kick 300 km s ^ 


Same as Table 2 but for kick 300 km s ^ model. 



underlOO 

overlOO 

140 

NS-NS 

8 (112) 

8 (112) 

4(78) 

NS-BH 

11922 (13133) 

11941 (13152) 

12115 (13330) 

BH-BH 

70728 (75011) 

78058 (82496) 

86876 (93481) 

merging NS-NS 

1 (85) 

1 (85) 

1 (60) 

merging NS-BH 

3893 (4483) 

3900 (4490) 

4406 (5002) 

merging BH-BH 

51928 (56021) 

58793 (63041) 

64084 (70252) 


Table 9. oA = 0.01 


Same as Table 2 but for a\ = 0.01 model. 



underlOO 

overlOO 

140 

NS-NS 

0 (0) 

0 (0) 

0 (0) 

NS-BH 

148290 (148770) 

148290 (148770) 

116548 (117117) 

BH-BH 

340893 (352047) 

345140 (363191) 

365526 (382686) 

merging NS-NS 

0 (0) 

0 (0) 

0 (0) 

merging NS-BH 

0 (294) 

0 (294) 

30 (412) 

merging BH-BH 

32283 (43437) 

36530 (54581) 

27790 (44950) 
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Table 10. a\ = 0.1 


Same as Table 2 but for a\ = 0.1 model. 



underlOO 

overlOO 

140 

NS-NS 

0 (0) 

0 (0) 

0 (0) 

NS-BH 

162814 (173016) 

162814 (173016) 

130556 (138835) 

BH-BH 

434590 (464369) 

448847 (480217) 

480520 (520031) 

merging NS-NS 

0 (0) 

0 (0) 

0 (0) 

merging NS-BH 

45 (181) 

45 (181) 

1065 (1877) 

merging BH-BH 

111696 (141356) 

125953 (157204) 

124830 (164240) 


Table 11. oA = 10 


Same as Table 2 but for aX = 10 model. 


underlOO 

overlOO 

140 

NS-NS 

NS-BH 

BH-BH 

merging NS-NS 
merging NS-BH 
merging BH-BH 

1116 (2215) 
198408 (198758) 
542399 (542603) 
890 (1949) 

767 (975) 
91787 (91989) 

1116 (2215) 
198408 (198758) 
560156 (560360) 
890 (1949) 

767 (975) 
104656 (104858) 

840 (1616) 
166173 (166408) 
624631 (624958) 
634 (1381) 

506 (645) 

93729 (94055) 


Table 12. /3 = 0 


Same as Table 2 but for /3 = 0 model. 



underlOO 

overlOO 

140 

NS-NS 

0 (279) 

0 (279) 

0 (195) 

NS-BH 

185335 (187638) 

185335 (187638) 

153435 (155694) 

BH-BH 

517067 (522581) 

534693 (540316) 

595894 (604930) 

merging NS-NS 

0 (279) 

0 (279) 

0 (195) 

merging NS-BH 

50 (149) 

50 (149) 

825 (1255) 

merging BH-BH 

115056 (120532) 

131060 (136645) 

128894 (137903) 


Table 13. /3 = 0.5 


Same as Table 2 but for /? = 0.5 model. 



underlOO 

overlOO 

140 

NS-NS 

5 (380) 

5 (380) 

6 (272) 

NS-BH 

193921 (196094) 

193921 (196094) 

158518 (160442) 

BH-BH 

549893 (554150) 

554966 (559228) 

628253 (635698) 

merging NS-NS 

5 (380) 

5 (380) 

6 (272) 

merging NS-BH 

199 (286) 

199 (286) 

766 (1082) 

merging BH-BH 

117094 (121310) 

119758 (123979) 

126090 (133512) 
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Table 14. /3 = 1 


Same as Table 2 but for /? = 1 model. 


underlOO overlOO 140 


NS-NS 

NS-BH 

BH-BH 

merging NS-NS 
merging NS-BH 
merging BH-BH 


1359 (2006) 
218311 (220521) 
531452 (536579) 
1358 (2005) 
119 (255) 
50119 (55214) 


1359 (2006) 
218311 (220522) 
531484 (536611) 
1358 (2005) 
119 (255) 
50119 (55214) 


898 (1344) 
178444 (180375) 
610732 (619230) 
898 (1344) 
578 (917) 
57025 (65121) 


Table 15. Worst 


Same as Table 2 but for Worst model. 



underlOO 

overlOO 

140 

NS-NS 

1637 (1637) 

1637 (1637) 

1604 (1604) 

NS-BH 

4345 (4345) 

4345 (4345 

4283 (4285) 

BH-BH 

5227 (5235) 

5227 (5235) 

5560 (5586) 

merging NS-NS 

1562 (1562) 

1562 (1562) 

1532 (1532) 

merging NS-BH 

1645 (1645) 

1645 (1645) 

1604 (1606) 

merging BH-BH 

3195 (3203) 

3195 (3203) 

3376 (3399) 
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Figure 1. our standard model 

Each line is the normalized distribution of the BH-BH chirp mass. 
The red, green, blue, pink, light blue and grey lines are the un¬ 
derlOO case with optimistic core-merger criterion, the overlOO 
case with optimistic core-merger criterion, the 140 case with opti¬ 
mistic core-merger criterion, the underlOO case with conservative 
core-merger criterion, the overlOO case with conservative core¬ 
merger criterion and the 140 case with conservative core-merger 
criterion, respectively. Ntotal = 10^ binaries. 
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Figure 2. IMFdogflat 
Same as Fig.l but for IMFdogfiat model. 
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Figure 3. IMF:Salpeter 
as Fig.l but for IMF:Salpeter model. 
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Figure 4. lEFxonst. 
as Fig.l but for IEF:const. model. 
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Figure 5. IEF:e 
as Fig.l but for lEFie"*^'^ model. 
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Figure 6. kick 100 km s ^ 
Same as Fig.l but for kick 100 km model. 


Figure 9. aX = 0.1 
Same as Fig.l but for aX = 0.1 model. 



20 40 60 80 100 120 140 

BH-BH chirp mass 



20 40 60 80 100 120 140 

BH-BH chirp mass [Mg^jJ 


Figure 7. kick 300 km s ^ 
Same as Fig.l but for kick 300 km s~^ model. 


Figure 10. aX = 10 
Same as Fig.l but for aX = 10 model. 
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Figure 8. aX = 0.01 


Figure 11. /3 = 0 


Same as Fig.l but for aX = 0.01 model. 


Same as Fig.l but for /3 = 0 model. 
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Same 


chirp mass [Mg, 
Figure 12. /3 = 0.5 
as Fig.l but for /3 = 0.5 model. 
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Figure 13. /3 = 1 

as Fig.l but for /3 = 1 model. 
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Figure 14. Worst 
as Fig.l but for Worst model. 



Figure 15. merger timeilMF 

This figure describes the merger time distributions of Pop III 
BH-BHs. The red line, the green line and the blue line are 
our standard model, the logfiat model and the Salpeter model. 
Ntotal = 10® binaries. 
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logt 

Figure 16. merger time:e 

This figure describes the merger time distributions of Pop III BH- 
BHs. The red line, the pink line and the light blue line are our 
standard model, IEF:const. model and IEF:e“^'^ model. Ntotal = 
10® binaries. 
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Figure 17. merger timeikick 

This figure describes the merger time distributions of Pop III 
BH-BHs. The red line, the orange line and the black line are our 
standard model, the kick 100 km model and the kick 300 km 
model. Ntotal = 10^ binaries. 



Figure 19. merger time:/3 

This figure describes the merger time distributions of Pop III BH- 
BHs. The red line, the green line, the blue line and the pink line 
are our standard model, the /3 = 0 model, the (3 = 0.5 and the 
jd = 1 model. Ntotal = 10^ binaries. 



Figure 18. merger time:aA 

This figure describes the merger time distributions of Pop HI BH- 
BHs. The red line, the orange line, the grey line and the black 
line are our standard model, the olX = 0.01 model, the qA = 0.1 
and the aX = 10 model. Ntotal = 10^ binaries. 


3.2 The merger rate properties of Pop III 
BH-BHs 


In this sub-section, we show the merger rate density of each 
model. The merger rate density Ri{t) is calculated using 
results of each model and the Pop III star formation rate of 
de Souza, Yoshida & loka (2011) (See Fig. 21 and Kinugawa 
et al. (2014) for details) as 


Riit) = 


It 


fb SFR{t') Ni{t - t') 


dt' 


(19) 


+ tb < M > A^total 
where i is the type of compact binaries such as NS-NS, NS- 
BH or BH-BH. /(, is the initial binary faction. The recent 
cosmological hydrodynamics simulation (Susa et al. 2014) 
suggests that the binary fraction is about 50%. Thus, we 
use fb = 1/2 since the total number of the binary is half of 
the total number of the stars in the binary. < M > is the 
mean initial stellar mass of Pop III star derived from IMF 



logt 

Figure 20. merger time:Worst 

This figure describes the merger time distributions of Pop III BH- 
BHs. The red line, the blue line, the light blue line, the black line, 
the orange line, the pink line and the green line are our standard 
model, the Salpeter model, the IEF:e“'^ ® model, the a\ = 0.01 
model, the /3 = 1 model and the Worst model. Ntotal = 10® 
binaries. 

and the initial mass range. SFR(T) [Mq yr“^ Mpc“®] is the 
Pop III star formation rate at t'. Ni{t — t') is the number of 
type i compact binaries which are formed in \t' + dt] and 

merge at time t. A/otai is the total number of the simulated 
binaries. 

Fig. 22 and not show the merger rate densities 
[Myr“^ Mpc“®] of BH-BHs as a function of cosmic time 
(lower abscissa) and redshift 2 (upper abscissa) in our stan¬ 
dard model and the worst model. It is seen that in each 
model the merger rate densities for the same redshift de¬ 
pend on neither the initial mass range ( [10,100] or [10,140] 
) nor the CE merger criterion. The other models have the 
same dependencies so that we do not show their figures. 
As a function of the redshift, the merger rate densities are 
nearly constant from 2 = 0 to 2 ~ 1 in each model. Ta- 
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z 

Figure 21. The star formation rate density (comoving) calculated by de Souza et al. (2011). The unit of the rate is Mq per comoving 
volume per proper time. The red line is the the total SFR density of Pop III stars. 


ble 16 shows the merger rate density [Myr”*^ Mpc“®] at 
z = 0 (tnubbie = 13.8 Gyrs) for each model. The lowest rate 
is as expected in worst model while the highest rate is in 
P = 0.5 model. 

Fig. 24 shows the difference of the merger rate density 
of each model for underlOO case. Table 17 describes the peak 
redshift of the BH-BHs merger rate density of each model 
in underlOO case. It is seen that the peak redshift of the 
BH-BHs merger rate density ranges from 8.8 to 7.15. These 
peak redshifts are near the peak of the star formation rate 
at z ~ 9. In the following, we discuss the difference of each 
model. 

The IMF dependence of the peak redshift of the merger 
rate density is clear seen. Namely for the steeper IMF the 
peak redshift is small although the difference is not so large 
(^0.45 in z). Since BH-BH progenitors whose initial mass is 
lower than 50 Mq tend to evolve via the RLOF but not via 
the CE, the steeper IMF can make BH-BH progenitors to 
evolve via RLOFs. BH-BHs which evolved via RLOFs tend 
to have the wider orbit than BH-BHs which evolved via CE 
phases. Therefore, the typical merger time for the steeper 
IMF tends to be long so that the peak redshift is smaller. 
As for the lEF dependence, no tendency is seen while as for 
the natal kick velocity dependence, the peak redshift for 100 
km s”*^ model is smaller than our standard model, but that 
of the 300 km s“^ model is higher than our standard model. 
In the 100 km s“^ model, the kick makes the BH-BHs to ec¬ 
centric orbit so that the merger time becomes smaller than 
that of the circular orbit. Thus, the number of the merging 
BH-BHs which merge at the high redshift should increase. 
However, the BH-BHs which cannot merge in our standard 
case due to wide orbit can merge due to the natal kick. Con¬ 
sequently, the number of the merging BH-BHs which merge 



The red, green, blue, pink, light blue and grey lines are the un¬ 
derlOO case with optimistic core-merger criterion, the overlOO 
case with optimistic core-merger criterion, the 140 case with opti¬ 
mistic core-merger criterion, the underlOO case with conservative 
core-merger criterion, the overlOO case with conservative core¬ 
merger criterion and the 140 case with conservative core-merger 
criterion, respectively. 

at lower redshift tends to increase. In the 300 km s“^ model, 
the natal kick velocity is too large so that the binary tends 
to disrupt. However, if the natal kick direction is against the 
orbital direction, the natal kick behaves like the brake of the 
car or if the separation before the natal kick is very close, 
the binary can survive. Thus, the survived binary tends to 
have very close and eccentric orbit so that they can merge 
early. This explains the apparent strange behavior of the 
dependence of the peak redshift on the natal kick velocity. 
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Table 16. The merger rate density [Myr ^ Mpc at 2 : = 0(tHubble = 13-8 Gyrs) 


underlOO overlOO 140 


our standard 
IMFdogflat 
IMFiSalpeter 
lEFxonst. 
IEF:e-0S 
kick 100 km s“^ 
kick 300 km s“^ 
aA = 0.01 
a\ = 0.1 
a\ = 10 

/3 = 0 
jS = 0.5 
/3 = 1 
Worst 


0.0258 (0.0260) 
0.0230 (0.0232) 
0.0116 (0.0117) 
0.0267 (0.0267) 
0.0252 (0.0252) 
0.0210 (0.0212) 
0.00726 (0.00732) 
0.00542 (0.00542) 
0.0249 (0.0255) 
0.0229 (0.0229) 
0.0258 (0.0260) 
0.0361 (0.0362) 
0.0186 (0.0187) 
0.00194 (0.00194) 


0.0277 (0.0279) 
0.0240 (0.0245) 
0.0121 ( 0 . 0122 ) 
0.0288 (0.0288) 
0.0270 (0.0271) 
0.0223 (0.0226) 
0.00747 (0.00754) 
0.00542 (0.00542) 
0.0249 (0.0255) 
0.0253 (0.0253) 
0.0277 (0.0279) 
0.0369 (0.0370) 
0.0186 (0.0187) 
0.00194 (0.00194) 


0.0251 (0.0252) 
0.0232 (0.0236) 
0.0131 (0.0132) 
0.0242 (0.0242) 
0.0228 (0.0228) 
0.0203 (0.0207) 
0.00657 (0.00672) 
0.00290 (0.00290) 
0.0207 (0.0210) 
0.0192 (0.0192) 
0.0251 (0.0252) 
0.0320 (0.0321) 
0.0159 (0.0161) 
0.00169 (0.00169) 
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Figure 23. Worst 
The same as Fig. 22 but for Worst model. 

In the case of the CE parameter, it changes the number 
of survived binaries during the CE phase and the merger 
time of the BH-BHs. If aA is low, i.e. the orbit energy loss 
via the CE phase is high, the number of survived binaries 
during the CE phase is small and the merger time of the 
BH-BHs is short. In qA = 0.01 model, the parameter aA 
is so small that the almost all binaries which enter the CE 
phase merge during the CE phase. Thus, merging BH-BH 
progenitors evolved via RLOF so that they have wide orbit. 
Therefore, their merger time tend to be long and the peak 
redshift is smaller than our standard model. In aA = 0.1 
model, aA is small but the binaries which become the CE 
phase do not merge during the CE phase and they have close 
orbit due to small aA. Thus, the merger time of BH-BHs is 
short and the peak redshift is large. In aA = 10 model, aA 
is so large that binaries after the CE phase have wide orbit 
due to large aA. Thus, the merger time of BH-BHs is long 
and the peak redshift is small. These consideration explains 
the strange behavior of the peak redshift on aA parameter. 

In the case of the parameter (5, it not only changes the 
mass accretion to the secondary but also changes the crite- 


Table 17. The peak redshift of the BH-BHs merger rate density 


model 

peak redshift 

our standard 

7.85 

IMFdogflat 

7.75 

IMFiSalpeter 

7.4 

lEFxonst. 

7.85 

IEF:e-'3® 

7.8 

kick 100 km 8“^ 

7.5 

kick 300 km 

8.65 

aA = 0.01 

7.2 

qA = 0.1 

8.5 

aA = 10 

6.85 

/3 = 0 

7.85 

P = 0.5 

7.15 

P = 1 

00 

bo 

Worst 

8.3 


rion of the dynamically unstable mass transfer. In /3 = 0.5 
model, the mass transfer is dynamically stable so that the 
number of the binaries which evolve via RLOF but not via 
the CE phase. Thus, the typical merger time is long and the 
peak redshift is low. In /I = 1 model, the mass transfer is 
always dynamically stable. Furthermore, the mass accretion 
to the secondary during RLOF does not occur so that the 
orbit becomes wide, On the other hand, binaries which have 
the eccentric orbit have only the CE phase. Thus, the merg¬ 
ing BH-BHs are separated into two groups. In one group 
the binaries evolve via RLOF while in the other group they 
evolve via the CE phase due to the eccentric orbit. The for¬ 
mer group merges at low redshift and the latter group does 
at high redshift. Therefore, in the /3 = 1 model, the merger 
rate density are bimodal as shown in Fig. 24. 

Note that the maximum merger rate density of 
3.7 X 10“® events yr“^ Mpc“® from Table 16 with 
P = 0.5 and overlOO case is consistent with the up¬ 
per limit of ~ 10“^ events yr“^ Mpc“® by LIGO- 

Virgo(S6/VSR2/VSR3) (Aasi et al. 2013). 


Z 

54 3 2 1 0.5 0.1 0 

underlOO (optimistic core merger) 
overlOO (optimistic core merger) ■■■■' 

140 (optimistic core merger) . 

underlOO (conservative core merger). . 

overlOO (conservative core merger) 

140 (conservative core merger). 
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Figure 24. The merger rate densities 


3.3 The detection rate of Pop III BH-BHs by the 
second generation detectors 

In this section, we show how to calculate the detection rate 
of BH-BHs. Our Pop HI population synthesis simulations 
produced a set of merging BH-BHs with component masses 
and merger time. To estimate the number of the event and 
the parameter decision accuracy according to the generated 
binary mass distribution, distances, various incident angles 
and orientations of the orbit plane, we employ the simple 
Monte-Carlo simulation. As the typical 2nd generation de¬ 
tectors that has a little advantage in a low frequency band in 
underground site, we employ KAGRA detector in our sim¬ 
ulation. We use the detection range of Kanda & the LCGT 
collaboration (2011). The official sensitivity limit of the KA- 
GRA^ is suitable for the detection of both inspiral and ring- 
down gravitational waves from the 10-30 Mq binaries. In 
the Monte-Carlo simulation, we placed each event at ran¬ 
dom position in the hemisphere, random direction of the bi¬ 
nary orbit plane. The direction in the cosmological redshift 
and the mass are given by the Pop HI binary simulation. 
We iterate many events for 1000 years, then, we estimate 
the expected detection rate for one year observation. The 
error of the rate is given as the square root of the number 
of detection in Monte-Carlo trials. 

Tables 18 to 23 describe the detection rate of BH-BHs 
in each model. This table shows the detection rates of Pop 
HI BH-BHs for underlOO cases with the optimistic core¬ 
merger criterion. The first column shows the name of the 
model. The second column shows the detection rate only 
by the inspiral chirp signal. The third, the fourth and the 
fifth columns show the detection rate only by the quasi nor¬ 
mal mode (QNM) with Kerr parameter a/M = 0.70, the 
detection rate by the quadrature sum of the inspiral chirp 
signal and the QNM with a/M = 0.70 and the detection 
rate by the linear sum of the inspiral chirp signal and the 
QNM with a/M = 0.70, respectively. The sixth, the sev¬ 
enth and the eighth columns show the detection rates only 
by the QNM with a/M = 0.98, the detection rate by the 

^ http://gwcenter.icrr.u-tokyo.ac.jp/en/researcher/parameter 


quadrature sum of the inspiral chirp signal and the QNM 
with a/M = 0.98 and the detection rate by the linear sum 
of the inspiral chirp signal and the QNM with a/M = 0.98, 
respectively. When signal-to-noise ratio of event that is cal¬ 
culated by matched filtering equation, over threshold S/N 
= 8, the event is detected. The QNM S/N is calculated by 
Eq. B14 in Flanagan & Hughes (1998). tr in this equation 
is the fraction of binary total mass energy radiated in the 
QNM. We assumed the value er = 0.03. Since their equation 
is averaged over the GW polarization and the sky location, 
a factor 1/5 is multiplied by the equation. However we have 
to take account of angular values of binary, we replace the 
factor with y^(l -|- cos^ '•)^/4 • + cos^ where the t 

is the inclination angle, and F+,Fx are KAGRA antenna 
pattern functions. For the fourth and sixth columns, their 
S/N are calculated by the linear summation of S/N of the 
inspiral and the QNM with a/M = 0.70 and 0.98, respec¬ 
tively. All the rates are based on 1000 years Monte Carlo 
simulations. 

In Zlochower & Lousto (2015), the reasonable value 
of Kerr parameter a/M is about 0.7. Thus, we focus on 
the detection rate by the quadrature sum of the inspiral 
chirp signal and the QNM with a/M = 0.70. The rates 
of the quadrature sum of the inspiral and the QNM is 
about 1/2 of the rates of the linear sum of the inspiral and 
the QNM. In our standard model with underlOO case and 
optimistic core-merger criterion, the detection rate by the 
quadrature sum of the inspiral chirp signal and the QNM is 
~1.8 xlO^ events yr“^ (SFRp/(10“^ ® Mq yr“^ Mpc“®)) • 
([fb/(l -I- fb)]/0.33) where SFRp and fb are the peak value 
of the Pop HI star formation rate and the binary fraction. 
Errsya = 1 corresponds to the rate for our standard model 
with underlOO case and the optimistic core-merger criterion. 
The definition of Errgys is slightly different from that in our 
previous paper (Kinugawa et al. 2014). That is, the new def¬ 
inition is based on Monte Carlo simulations of the detection 
of the inspiral chirp signal and QNM with a/M = 0.7. The 
basic numerical data of population synthesis of our standard 
model is the same. Then Tables 18-23 show that Errsys 
ranges from 4.6 x 10“^ to 4 for a/M = 0.7. This means that 
the detection rate of the coalescing Pop HI BH-BHs ranges 
8.3 —7.2x 10^ events yr“^ (SFRp/(10“^ ®Mq yr"*^ Mpc“®)) 
■([fb/(l-f fb)]/0.33). The minimum detection rate of 
the coalescing Pop HI BH-BHs corresponds to the 
worst model which we think unlikely so that unless 
(SFRp/(10-" ® Mq yr-i Mpc'^)) • ([fb/(l -f fb)]/0.33) < 0.1, 
we can expect the Pop HI BH binary merger at least one 
event per year by the second generation gravitational wave 
detector. 

Fig. 3 of Nakano et al. (2015) shows that if the S/N 
of the Pop HI BH-BH QNM is 50, we can check the gen¬ 
eral relativity with the significance of much more than 5 
sigma level. The criterion of S/N whether the significance 
has more than 5 sigma is 35. Therefore, we expect the 
enough accuracy to discuss about the GR test with at least 
one event of S/N= 35. In our standard model, the detec¬ 
tion rate of Pop HI BH-BHs whose S/N is more than 35 is 
3.2 events yr-i(SFRp/(10-^ ® Mq yr’^ Mpc”®)) • ([fb/(l -f 
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fb)]/0.33) • Errsys- Thus, there might be a good chance to 
check whether the GR is correct or not in the strong gravity 
region. The forth and the seventh columns of Tables 24-29 
show that there is a good chance to observe QNM of the 
merged BH since it reflects the space-time near the BH and 
their complex frequency does not depend on how it is ex¬ 
cited. The detection of the expected QNM of the BH can 
confirm the GR in the strong gravity region. If it is differ¬ 
ent from the expected value, the true theory of gravity is 
different from GR. 


4 DISCUSSION & SUMMARY 

In this paper, we performed the Pop III binary population 
synthesis and examined the parameter dependence of Pop 
III binary evolutions. We examined the dependence of the 
results on IMF, lEF, the natal kick velocity, the CE pa¬ 
rameters and the lose fraction of stellar mass. As for the 
chirp mass distribution, each model has the peak at around 
30 Mq. In several models, the chirp mass distribution has a 
tail from 30 Mq to more massive region. However the robust 
property is that the chirp mass distribution has the peak at 
30 Mq. 

In order to compare the variability of Errsys, we refer 
previous researches such as Belczynski et al. (2012); Do- 
minik et al. (2015). In Belczynski et al. (2012), they cal¬ 
culated the solar metal [Z = 0 . 02 ) binaries and 10 % solar 
metal {Z = 0.002) binaries to estimate the detection rates 
assuming as half of the stars formed with solar metal and the 
other withl0% solar metal. In Belczynski et al. (2012), they 
calculated 20 models by varying the maximum NS mass, 
the natal kick velocity, rapid or delayed supernova models 
which change the mass spectrum of supernova remnants, 
wind mass loss and /3. They also considered whether in the 
Hertzsprung gap donors always merge companion during the 
CE phase or not. The detection rate of their realistic Stan¬ 
dard model in Belczynski et al. (2012) is 517.3 events yr~^ 
The detection rates of Belczynski et al. (2012) are from 
14 events yr“^ to 12434.4 events yr“^. Thus, Errsys of Bel¬ 
czynski et al. (2012) is from 2.7 x 10“^ to 24. 

On the other hand, Dominik et al. (2015) calculated 
binaries whose metallicity range is from Z = 10“^ to 
Z = 0.03 and estimated the detection rate using the metal¬ 
licity and SFR evolution models. There are 16 models vary¬ 
ing high-end or low-end metallicity models, whether in the 
Hertzsprung gap donors always merge companion during 
the CE phase or not, rapid or delayed supernova models, 
the natal kick and waveform models. The detection rate 
of their Standard model of high-end metallicity scenario 
in Dominik et al. (2015) is 306 events yr“^ The detection 
rate of Dominik et al. (2015) is from 8.2 events yr“^ to 
3087 events yr“^ by the 3-detector network using inspiral 
and PhC waveform (S/N=10). Thus, Errsys of Dominik et 
al. (2015) is from 2.7 x 10“^ to 10. Therefore, the variabil¬ 
ity of Errsys of Pop HI is less than that of Pop I and Pop 
H, although the models are different. There are two reasons 
for this difference. Firstly, Pop HI star binaries do not en¬ 
ter the CE phase so that the result does not depend on the 


treatment of the CE phase so much. Secondly, the Pop HI 
compact binaries are more massive than the Pop I compact 
binary so that the Pop HI binaries are hard to be disrupted 
by the natal kick. Therefore, the property of the chirp mass 
distribution and the detection rate are likely robust result. 
However, note that in the case of the detection rate of Pop 
HI there are the dependence on the SFR and fb yet. 

There are some uncertainties yet such as the separa¬ 
tion distribution function and the mass ratio distribution 
function for which we did not alter. The former will change 
the number of close binaries which can have binary interac¬ 
tions. Therefore this effect may change the event rate, but 
the property of chirp mass distribution is not likely changed 
a lot because the binary interaction is not changed. From 
our Monte Cairo simulations, the chirp mass distribution of 
Pop HI BH-BHs is upward to the high mass and has a peak 
at ~ 30 Mq in each model. The compact objects in ICIO 
X-1 and NGC300 X-1 may be around 30 M© and they might 
become coalescing massive BH-BHs whose chirp masses are 
11-26 M 0 (See Bulik et al. (2011)). Thus, Pop I stars or Pop 
H stars might become coalescing massive BH-BHs. However, 
the observed typical mass of Pop I BH-BHs is around 10 M© 
and massive BH like ICIO X-1 and NGC300 X-1 would be 
rare (See also Fig. 1 in Belczynski et al. (2012)) so that the 
chirp mass distribution of Pop I BH-BHs might be flat or 
decreasing as a function of mass. The result of the binary 
population synthesis simulation for Pop I and Pop H stars 
by Dominik et al. (2015) also suggests that the chirp mass 
distribution of Pop I BH-BHs might be flat or downward to 
high mass (See the Fig. 7 in Dominik et al. (2015)). Fur¬ 
thermore, the Pop I and Pop H BH-BH detection rate of 
the standard model in Dominik et al. (2015) is 306 yr“^. 
The fraction of the Pop I and Pop II BH-BH whose mass is 
larger than 20 Mq is about 25%. Thus, the detection rate 
of the Pop I and Pop H high mass BH-BHs is expected as 
about 80 yr“^. Note that this value depends on the Errsys 
of Dominik et al. (2015) which is from 2.7 x 10~^ to 10. 
Therefore, if the detection rate of the coalescing Pop I and 
Pop H high mass BH-BHs is lower than that of Pop HI, 
we may be able to confirm the existence of Pop HI star by 
the detection of the chirp signal and QNM to determine the 
chirp mass and the total mass distribution since the typical 
mass of Pop HI BH-BH binary is much larger than those 
of observed Pop I BH. On the other hand, if the detection 
rate of the coalescing Pop I and Pop H high mass BH-BHs 
is higher, Pop HI BH-BHs contribute only some parts of 
the gravitational wave events of BH-BHs. In this case, the 
existence of Pop HI binaries will be conhrmed by the in¬ 
vestigation of the merger rate history as function of redshift 
by DECIGO (DECi hertz Interferometer Gravitational wave 
Observatory) (Seto, Kawamura & Nakamura 2001). 

As for the mass ratio distribution function, if the num¬ 
ber of the high mass ratio (i.e. near 1 ) increases, the num¬ 
ber of BH-BH probably increase. On the other hand, if the 
number of the low mass ratio increase, the number of BH- 
BH will decrease while the number of NS-BH will increase. 
We will check the dependence of these two initial distribu¬ 
tion functions in future work. Development in the simulation 
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Table 18. underlOO cases with optimistic core-merger criterion, 1000 years, S/N ^ 8 

This table shows the detection rates of Pop III BH-BHs for underlOO cases with the optimistic core-merger criterion. The first column 
shows the name of the model. The second column shows the detection rate only by the inspiral chirp signal. The third, the fourth and 
the fifth columns show the detection rate only by the quasi normal mode (QNM) with Kerr parameter a/M = 0.70, the detection rate 
by the quadrature sum of the inspiral chirp signal and the QNM with a/M = 0.70 and the detection rate by the linear sum of the 
inspiral chirp signal and the QNM with a/M = 0.70, respectively. The sixth, the seventh and the eighth columns show the detection 
rates only by the QNM with a/M = 0.98, the detection rate by the quadrature sum of the inspiral chirp signal and the QNM with a/M 
= 0.98 and the detection rate by the linear sum of the inspiral chirp signal and the QNM with a/M = 0.98, respectively. When 
signal-to-noise ratio of event that is calculated by matched filtering equation, over threshold S/N = 8, the event is detected. All the 

rates are based on 1000 years Monte Carlo simulations. 


14models 

Inspiral 
S/N^ 8 
[/lOOO yrs] 

QNM(0.70) 
S/N^ 8 
[/lOOO yrs] 

Quadrature sum of 
Inspiral and 
QNM(0.70) 
S/Nls 8 
i/iooo yrs] 

Linear sum of 
Inspiral and 
QNM(0.70) 
S/N^ 8 
1/1000 yrs] 

QNM(0.98) 
S/Nls 8 
1/1000 yrs] 

Quadrature sum of 
Inspiral and 
QNM(0.98) 
S/N5: 8 

1/1000 yrs] 

Linear sum of 
Inspiral and 
QNM(0.98) 
S/N^ 8 
1/1000 yrs] 

our standard 

85747 

67337 

180086 

392542 

10680 

110896 

213115 

IMFdogflat 

74764 

43130 

139705 

305428 

6524 

91378 

167958 

IMF:Salpeter 

47055 

16153 

74809 

156804 

2201 

53960 

90954 

lEFxonst. 

80947 

65225 

173059 

380515 

10476 

105409 

205489 

IEF:e-°® 

77922 

63050 

167289 

369763 

10342 

101941 

198808 

kick 100 km 

66901 

57590 

148370 

330875 

9243 

88399 

176161 

kick 300 km 8“^ 

16962 

14166 

36655 

80711 

2303 

22219 

43214 

qA=0.01 

18900 

13327 

37755 

82585 

2160 

23802 

44548 

«A=0.1 

67189 

73207 

165215 

366096 

11375 

93172 

193767 

aA=10 

74276 

59901 

156419 

336621 

10585 

96410 

182207 

/3=0 

85578 

67608 

180362 

393241 

10990 

110901 

213343 

/3=0.5 

117738 

74535 

229389 

498512 

11015 

145968 

272362 

/3=1 

58823 

45926 

123468 

268790 

7353 

75646 

145760 

Worst 

4922 

2322 

8449 

17433 

316 

5786 

9997 


Table 19. overlOO 

cases with the optimistic core-merger criterion, 1000 years, S/N ^ 8 

The same as Table 18 but for overlOO cases with the optimistic 

core-merger 

criterion. 





Quadrature sum 

of Linear sum of 


Quadrature sum of 

Linear sum of 




Inspiral and 

Inspiral and 


Inspiral and 

Inspiral and 


Inspiral 

QNM(0.70) 

QNM(0.70) 

QNM(0.70) 

QNM(0.98) 

QNM(0.98) 

QNM(0.98) 


S/N^ 8 

S/N^ 8 

S/Nls 8 

S/NIJ 8 

S/Nls 8 

S/N^ 8 

S/N^ 8 

14models 

1/1000 yrs] 

[/lOOO yrs] 

1/1000 yrs] 

1/1000 yrs] 

1/1000 yrs] 

1/1000 yrs] 

1/1000 yrs] 

our standard 

84174 

218283 

332447 

611583 

41503 

142496 

275799 

IMFdogflat 

74305 

130893 

229061 

435444 

24041 

109695 

204510 

IMF:Salpeter 

47216 

39218 

98543 

191702 

6890 

59307 

101093 

lEFxonst. 

79623 

222826 

332472 

609477 

42784 

139233 

271529 

IEF:e-0® 

77427 

221260 

327610 

598925 

43430 

136951 

266084 

kick 100 km s“l 

66563 

151555 

244415 

473699 

27246 

107999 

215263 

kick 300 km s“^ 

16843 

40830 

63663 

121315 

7334 

27728 

54613 

qA=0.01 

16864 

31946 

55038 

111547 

5170 

25307 

49459 

aA=0.1 

62295 

169260 

257276 

499826 

27164 

104855 

217745 

aA=10 

76876 

181513 

284367 

531845 

34200 

125362 

240192 

/3=0 

84596 

218993 

333565 

613502 

41842 

143224 

276753 

/3=0.5 

118459 

98205 

255012 

540234 

15457 

152208 

285242 

/3=1 

58887 

46219 

123729 

268745 

7482 

76071 

145729 

Worst 

4873 

2361 

8388 

17640 

330 

5698 

9974 
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Table 20. 140 cases with the optimistic core-merger criterion, 1000 years, S/N ^ 8 

The same as Table 18 but for 140 cases with the optimistic core-merger criterion 


14models 

Inspiral 

S/N5: 8 
[/lOOO yrs] 

QNM(0.70) 
S/N^ 8 
[/lOOO yrs] 

Quadrature sum of 
Inspiral and 
QNM(0.70) 
S/N> 8 
i/iooo yrs] 

Linear sum of 
Inspiral and 
QNM(0.70) 
S/n;s 8 

1/1000 yrs] 

QNM(0.98) 
S/N> 8 
1/1000 yrs] 

Quadrature sum of 
Inspiral and 
QNM(0.98) 
S/N5: 8 

1/1000 yrs] 

Linear sum of 
Inspiral and 
QNM(0.98) 
S/N^ 8 
1/1000 yrs] 

our standard 

57415 

526720 

605703 

880465 

229397 

301747 

468129 

IMFJogflat 

62822 

334314 

418581 

643524 

135840 

210816 

336500 

IMF:Salpeter 

45680 

94589 

153225 

262258 

34210 

85487 

139481 

lEFxonst. 

54288 

517651 

593218 

857973 

223233 

292210 

453214 

IEF:e-0® 

52191 

510059 

582984 

840165 

217287 

283573 

439967 

kick 100 km 

44081 

501665 

565394 

816311 

181426 

239150 

389811 

kick 300 km 

10582 

126851 

142451 

207343 

56333 

70251 

108314 

qA=0.01 

9265 

34743 

47728 

85457 

6163 

17465 

34075 

aA=0.1 

46225 

318804 

387640 

647921 

84572 

146699 

280102 

aA=10 

46837 

394562 

457263 

656125 

207658 

265181 

389217 

/3=0 

57037 

526812 

606147 

881259 

228443 

300762 

467266 

/9=0.5 

75890 

603167 

706657 

1032803 

262573 

355852 

552073 

13=1 

34158 

256049 

304003 

460519 

76025 

119493 

209309 

Worst 

4101 

4262 

9803 

20095 

701 

5460 

10240 


Table 21. underlOO 

cases with the conservative core-merger criterion, 1000 years, S/N ^ 8 

The same as Table 18 but for underlOO cases with the conservative core-merger 

criterion 





Quadrature sum of 

Linear sum of 


Quadrature sum of 

Linear sum of 




Inspiral and 

Inspiral and 


Inspiral and 

Inspiral and 


Inspiral 

QNM(0.70) 

QNM(0.70) 

QNM(0.70) 

QNM(0.98) 

QNM(0.98) 

QNM(0.98) 


S/N^ 8 

S/N^ 8 

S/N> 8 

S/n;s 8 

S/N> 8 

S/N^ 8 

S/N^ 8 

14models 

1/1000 yrs] 

[/lOOO yrs] 

[/lOOO yrs] 

[/lOOO yrs] 

[/lOOO yrs] 

[/lOOO yrs] 

[/lOOO yrs] 

our standard 

85078 

67694 

179779 

392884 

10940 

110427 

212302 

IMFJogflat 

75150 

43467 

140349 

306831 

6449 

91710 

168361 

IMF:Salpeter 

47239 

16174 

74706 

156639 

2235 

54032 

90505 

lEFxonst. 

81437 

65772 

173889 

382036 

10614 

106027 

206030 

IEF:e-0® 

77642 

63262 

167088 

369913 

10183 

101613 

198231 

kick 100 km s“l 

67689 

58360 

150101 

334191 

9400 

89472 

177738 

kick 300 km s“^ 

17426 

14650 

37703 

83042 

2355 

22851 

44532 

qA=0.01 

16187 

21897 

43744 

98315 

4138 

23912 

50245 

«A=0.1 

72430 

64944 

162394 

361198 

9863 

95990 

192033 

«A=10 

74185 

60266 

157326 

337812 

10294 

96638 

182698 

/3=0 

85590 

67951 

181242 

395187 

10850 

110854 

214213 

/3=0.5 

117864 

74766 

229642 

498925 

11113 

146208 

272465 

/3=1 

58961 

45682 

123120 

267733 

7462 

75987 

145467 

Worst 

4757 

2270 

8347 

17605 

298 

5603 

10014 


may make it possible to clarify initial conditions of Pop III 
binary. 

The Pop III star formation rate will determine the 
merger rate. Our using Pop III SFR (de Souza et al. 2011) is 
Fig. 21. There are some arguments on Pop III SFR besides 
de Souza et al. (2011). For example, Johnson et al. (2013) 
simulated the Pop III SFR by the smooth particle hydro¬ 
dynamics (SPH) simulations. In their simulations, the peak 
value of Pop III SFR is from ^ 10“® ’^ Mq yr“^ Mpc“® 


to ~ 10“® Mq yr“^ Mpc“® at 2 ~ 10. The difference be¬ 
tween the high value and the low value comes from without 
Lyman-Werner (LW) feedback or with LW feedback. Note 
that the result of these simulations might change if the metal 
pollution model changes. On the other hand, Kulkarni et al. 
(2014) and Yajima & Khochfar (2015) studied the Pop III 
SFR by considering the contribution of Pop III stars to cos¬ 
mic reionization. Kulkarni et al. (2014) suggest that the peak 
value of Pop III SFR is from 10“"^’^ Mq yr“^ Mpc“® to 
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Table 22. overlOO cases with the conservative core-merger criterion, 1000 years, S/N is 8 

The same as Table 18 but for overlOO cases with the conservative core-merger criterion 


14models 

Inspiral 

S/N5: 8 
[/lOOO yrs] 

QNM(0.70) 
S/N^ 8 
[/lOOO yrs] 

Quadrature sum of 
Inspiral and 
QNM(0.70) 
S/Nls 8 
i/iooo yrs] 

Linear sum of 
Inspiral and 
QNM(0.70) 
S/Nls 8 
1/1000 yrs] 

QNM(0.98) 
s/n;s 8 

1/1000 yrs] 

Quadrature sum of 
Inspiral and 
QNM(0.98) 
S/N5: 8 

1/1000 yrs] 

Linear sum of 
Inspiral and 
QNM(0.98) 
S/N^ 8 
1/1000 yrs] 

our standard 

84864 

214329 

329356 

608170 

41205 

142681 

275555 

IMFdogflat 

75062 

128962 

227785 

434489 

23497 

110212 

205048 

IMF:Salpeter 

47104 

38994 

98198 

191541 

6847 

59142 

101187 

lEFxonst. 

80521 

219273 

330324 

606041 

42209 

139338 

271312 

IEF:e-0® 

77318 

218702 

325574 

595975 

42196 

135838 

265493 

kick 100 km 

67736 

150146 

243859 

474159 

26286 

108305 

216333 

kick 300 km 

17139 

40582 

64028 

122219 

7137 

27864 

55062 

qA=0.01 

14194 

57924 

78827 

153302 

9617 

28113 

59507 

oA=0.1 

66817 

166466 

259438 

497156 

27926 

109303 

219794 

oA=10 

76476 

182640 

285088 

532376 

34371 

125343 

240240 

13=0 

84415 

213727 

328158 

607593 

40698 

141795 

274286 

13=0.5 

119110 

97920 

255288 

540990 

15356 

152569 

286062 

/3=1 

59051 

45600 

123150 

267169 

7314 

75779 

145191 

Worst 

4819 

2337 

8335 

17592 

333 

5702 

10007 


Table 23. 140 cases 

with the conservative core-merger criterion, 1000 years, S/N ^ 8 

The same as Table 18 but for 140 cases with the conservative 

core-merger 

criterion 





Quadrature sum of Linear sum of 


Quadrature sum of 

Linear sum of 




Inspiral and Inspiral and 


Inspiral and 

Inspiral and 


Inspiral 

QNM(0.70) 

QNM(0.70) QNM(0.70) 

QNM(0.98) 

QNM(0.98) 

QNM(0.98) 


S/N^ 8 

S/N^ 8 

S/Nis 8 S/Nls 8 

s/n;s 8 

S/N^ 8 

S/N^ 8 

14models 

1/1000 yrs] 

[/lOOO yrs] 

[/lOOO yrs] [/lOOO yrs] 

[/lOOO yrs] 

[/lOOO yrs] 

[/lOOO yrs] 

our standard 

55974 

538382 

616793 893170 

249666 

321038 

489786 

IMFdogflat 

62224 

337179 

420100 644264 

146827 

221150 

346650 

IMF:Salpeter 

45352 

97882 

155975 265370 

37519 

88623 

142980 

lEFxonst. 

53720 

528461 

603514 869645 

239380 

307730 

470677 

IEF:e-0® 

51664 

515375 

587375 845467 

229457 

295251 

452615 

kick 100 km s“l 

45815 

506994 

573221 830258 

184386 

244263 

399027 

kick 300 km s“^ 

11241 

135797 

152230 219985 

66444 

81421 

121657 

qA=0.01 

7431 

62666 

74189 123411 

14129 

24092 

47107 

aA=0.1 

47406 

366456 

437298 710775 

113247 

176176 

316298 

aA=10 

46402 

395992 

458184 656940 

210422 

267360 

391782 

(3=0 

55773 

539238 

617239 892843 

249140 

320150 

489898 

(3=0.5 

74932 

619682 

722387 1051940 

277705 

370639 

570146 

/3=1 

34014 

269616 

317328 480619 

79988 

123652 

216096 

Worst 

4276 

4663 

10282 20904 

981 

5888 

10703 


^ 10“^'® Mq yr“^ Mpc“® at 2 10. The difference be¬ 

tween the high value and the low value comes from that 
of the metal pollution timescale. While Yajima & Khoch- 
far (2015) suggests that the peak value of Pop III SFR is 
10“® Mq yr“^ Mpc“® at 2 15 in order to recover the 

observed Thomson scattering optical depth of the cosmic 
microwave background. The SFR of Pop III is controversial 
now. However, these estimated value of the SFR tell us that 


except for the worst model, we might expect the detection 
of GW from massive Pop III BH-BH near future. 

The present merger rate density which is calculated only 
from SFR between z=7 to z= 11 is about 50% of the whole 
merger rate density. Therefore SFRpeak is a good parameter 
in our adopted model of SFR. However, there is Pop HI SFR 
whose peak region is higher redshift than 2 10 such as the 

models in Yajima et al. (2015). In such a case, we have to 
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Table 24. underlOO cases with optimistic core-merger criterion, 1000 years, S/N ^ 35 

This table shows the detection rates of Pop III BH-BHs for underlOO cases with the optimistic core-merger criterion. The first column 
shows the name of the model. The second column shows the detection rate only by the inspiral chirp signal. The third, the fourth and 
the fifth columns show the detection rate only by the quasi normal mode (QNM) with Kerr parameter a/M = 0.70, the detection rate 
by the quadrature sum of the inspiral chirp signal and the QNM with a/M = 0.70 and the detection rate by the linear sum of the 
inspiral chirp signal and the QNM with a/M = 0.70, respectively. The sixth, the seventh and the eighth columns show the detection 
rates only by the QNM with a/M = 0.98, the detection rate by the quadrature sum of the inspiral chirp signal and the QNM with a/M 
= 0.98 and the detection rate by the linear sum of the inspiral chirp signal and the QNM with a/M = 0.98, respectively. When 
signal-to-noise ratio of event that is calculated by matched filtering equation, over threshold S/N = 35, the event is detected. All the 

rates are based on 1000 years Monte Carlo simulations. 


14models 

Inspiral 
S/N^ 35 
[/lOOO yrs] 

QNM(0.70) 
S/N5s 35 
[/lOOO yrs] 

Quadrature sum of 
Inspiral and 
QNM(0.70) 
S/N^ 35 
i/lOOO yrs] 

Linear sum of 
Inspiral and 
QNM(0.70) 
S/N> 35 
1/1000 yrs] 

QNM(0.98) 
S/N^ 35 
1/1000 yrs] 

Quadrature sum of 
Inspiral and 
QNM(0.98) 
S/N^ 35 
1/1000 yrs] 

Linear sum of 
Inspiral and 
QNM(0.98) 
S/N5: 35 
1/1000 yrs] 

our standard 

2166 

488 

3172 

7190 

66 

2395 

4052 

IMFdogflat 

1918 

326 

2591 

5544 

48 

2062 

3292 

IMF:Salpeter 

1087 

102 

1353 

2719 

13 

1150 

1711 

lEFxonst. 

2016 

476 

2945 

6679 

62 

2238 

3732 

IEF:e-°® 

1937 

433 

2855 

6362 

58 

2156 

3590 

kick 100 km 

1636 

404 

2456 

5623 

60 

1850 

3123 

kick 300 km 8“^ 

387 

91 

589 

1342 

9 

434 

732 

qA=0.01 

440 

96 

659 

1425 

14 

483 

832 

aA=0.1 

1795 

520 

2882 

6717 

69 

2050 

3569 

«A=10 

1829 

436 

2733 

5968 

72 

2027 

3407 

/3=0 

2118 

468 

3102 

7220 

59 

2345 

4011 

/3=0.5 

3176 

557 

4379 

9340 

52 

3447 

5512 

13=1 

1574 

352 

2315 

5057 

50 

1757 

2917 

Worst 

111 

21 

146 

297 

4 

118 

173 


Table 25. overlOO 

cases with the optimistic core-merger criterion, lOOOyears, S/N ^ 35 

The same as Table 24 but for overlOO cases with the optimistic 

core-merger 

criterion. 





Quadrature sum 

of Linear sum of 


Quadrature sum of 

Linear sum of 




Inspiral and 

Inspiral and 


Inspiral and 

Inspiral and 


Inspiral 

QNM(0.70) 

QNM(0.70) 

QNM(0.70) 

QNM(0.98) 

QNM(0.98) 

QNM(0.98) 


S/N^ 35 

S/N^ 35 

S/N> 35 

S/Nls 35 

S/Nls 35 

S/N^ 35 

S/N^ 35 

14models 

1/1000 yrs] 

1/1000 yrs] 

1/1000 yrs] 

1/1000 yrs] 

1/1000 yrs] 

1/1000 yrs] 

1/1000 yrs] 

our standard 

2122 

959 

3707 

8269 

244 

2587 

4693 

IMFdogflat 

1866 

543 

2867 

6241 

117 

2167 

3667 

IMF:Salpeter 

1098 

174 

1485 

2976 

36 

1202 

1881 

lEFxonst. 

1939 

953 

3475 

7756 

276 

2444 

4413 

IEF:e-0® 

1914 

928 

3464 

7738 

216 

2396 

4408 

kick 100 km 8“^ 

1677 

660 

2850 

6442 

175 

2008 

3642 

kick 300 km s“^ 

382 

172 

668 

1512 

40 

472 

863 

qA=0.01 

410 

158 

689 

1454 

35 

479 

839 

aA=0.1 

1692 

780 

3054 

6900 

159 

2066 

3923 

aA=10 

1762 

805 

3092 

6885 

197 

2171 

3900 

/3=0 

2247 

1002 

3808 

8397 

266 

2752 

4799 

/3=0.5 

2981 

637 

4349 

9394 

95 

3328 

5569 

[3=1 

1554 

341 

2294 

5071 

53 

1753 

2909 

Worst 

104 

13 

148 

331 

2 

116 

191 
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Table 26. 140 cases with the optimistic core-merger criterion, 1000 years, S/N ^ 35 

The same as Table 24 but for 140 cases with the optimistic core-merger criterion 


14models 

Inspiral 
S/N^ 35 
[/lOOO yrs] 

QNM(0.70) 
S/NS! 35 
[/lOOO yrs] 

Quadrature sum of 
Inspiral and 
QNM(0.70) 
S/Nls 35 
i/lOOO yrs] 

Linear sum of 
Inspiral and 
QNM(0.70) 
S/N> 35 
1/1000 yrs] 

QNM(0.98) 
S/Nls 35 
1/1000 yrs] 

Quadrature sum of 
Inspiral and 
QNM(0.98) 
S/NS! 35 
1/1000 yrs] 

Linear sum of 
Inspiral and 
QNM(0.98) 
S/NS! 35 
1/1000 yrs] 

our standard 

1410 

4487 

6553 

12149 

1017 

2757 

5342 

IMFdogflat 

1580 

2585 

4655 

8876 

572 

2375 

4371 

IMF:Salpeter 

955 

624 

1823 

3652 

137 

1185 

2018 

lEFxonst. 

1387 

4425 

6416 

11722 

954 

2677 

5252 

IEF:e-0® 

1318 

4220 

6120 

11353 

966 

2566 

5002 

kick 100 km 

1040 

2921 

4415 

8824 

735 

2028 

4072 

kick 300 km 

288 

968 

1370 

2489 

206 

538 

1064 

qA=0.01 

236 

139 

436 

933 

44 

318 

556 

oA=0.1 

1187 

1606 

3331 

7358 

397 

1823 

3789 

oA=10 

1210 

4366 

6072 

10415 

949 

2418 

4404 

/3=0 

1381 

4545 

6617 

12117 

962 

2718 

5292 

13 = 0.5 

1904 

4858 

7457 

14128 

956 

3223 

6230 

/3=1 

880 

1377 

2633 

5571 

393 

1470 

2877 

Worst 

86 

19 

115 

248 

3 

99 

147 


Table 27. underlOO 

cases with the conservative core-merger criterion, 1000 years, S/N ^ 35 

The same as Table 24 but for underlOO cases with the conservative core-merger 

criterion 





Quadrature sum of 

Linear sum of 


Quadrature sum of 

Linear sum of 




Inspiral and 

Inspiral and 


Inspiral and 

Inspiral and 


Inspiral 

QNM(0.70) 

QNM(0.70) 

QNM(0.70) 

QNM(0.98) 

QNM(0.98) 

QNM(0.98) 


S/NS! 35 

S/NS! 35 

S/N> 35 

S/Nls 35 

S/NlJ 35 

S/NS! 35 

S/NS! 35 

14models 

1/1000 yrs] 

1/1000 yrs] 

1/1000 yrs] 

1/1000 yrs] 

1/1000 yrs] 

1/1000 yrs] 

1/1000 yrs] 

our standard 

2109 

498 

3163 

7092 

54 

2350 

3986 

IMFdogflat 

1863 

302 

2646 

5530 

40 

2034 

3307 

IMF:Salpeter 

1132 

132 

1421 

2755 

17 

1200 

1804 

lEFxonst. 

2012 

493 

3016 

6707 

67 

2241 

3816 

IEF:e-0® 

1905 

438 

2792 

6365 

53 

2122 

3519 

kick 100 km s“l 

1641 

384 

2492 

5613 

61 

1825 

3131 

kick 300 km s“^ 

376 

80 

565 

1278 

11 

417 

720 

qA=0.01 

396 

164 

670 

1584 

16 

467 

833 

«A=0.1 

1962 

525 

2939 

6624 

69 

2191 

3714 

«A=10 

1757 

375 

2551 

5695 

52 

1913 

3193 

( 3=0 

2221 

521 

3213 

7150 

64 

2458 

4108 

( 3 = 0.5 

3168 

555 

4401 

9282 

73 

3454 

5532 

/3=1 

1506 

330 

2223 

4925 

48 

1675 

2843 

Worst 

80 

13 

123 

301 

1 

89 

170 


consider j SFR{z)dz in order to compare the dependence 
on the SFR. 

We discuss also the binary fraction of Pop III. Recently, 
the resolution of the multi-dimension simulation becomes so 
high that the fragmentation of disk at the Pop III stellar 
formation can be studied like in (Clark et al. 2011). The 
recent cosmological hydrodynamics simulation (Susa et al. 
2014) suggests that the binary fraction is about 50%. Thus, 
we use fb = 1/2 since the total number of the binary is 


half of the total number of the stars in the binary. However, 
the binary fraction is controversial. Thus, we express the 
uncertainty of the binary fraction of Pop III as fp. 

From Fig. 24, the peak of the rate of the merger of Pop 
III BH binary is around 2 ; = 9 so that the observed frequency 
of the chirp signal and the quasi normal mode are ~ 10 times 
small. To detect such a low frequency gravitational wave, 
DECIGO (Seto, Kawamura & Nakamura 2001) will be most 
appropriate. When DECIGO starts an observation around 
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Table 28. overlOO cases with the conservative core-merger criterion, 1000 years, S/N ^ 35 

The same as Table 24 but for overlOO cases with the conservative core-merger criterion 


14models 

Inspiral 
S/N^ 35 
[/lOOO yrs] 

QNM(0.70) 
S/N5: 35 
[/lOOO yrs] 

Quadrature sum of 
Inspiral and 
QNM(0.70) 

S/N> 35 
i/lOOO yrs] 

Linear sum of 
Inspiral and 
QNM(0.70) 
S/N> 35 
1/1000 yrs] 

QNM(0.98) 
S/Nis 35 
1/1000 yrs] 

Quadrature sum of 
Inspiral and 
QNM(0.98) 

S/N^ 35 

1/1000 yrs] 

Linear sum of 
Inspiral and 
QNM(0.98) 
S/N5: 35 
1/1000 yrs] 

our standard 

2103 

936 

3682 

8228 

243 

2543 

4675 

IMFdogflat 

1918 

598 

2945 

6221 

155 

2231 

3697 

IMF:Salpeter 

1107 

210 

1483 

2977 

42 

1215 

1879 

lEFxonst. 

1970 

884 

3410 

7700 

224 

2436 

4301 

IEF:e-0® 

1881 

922 

3364 

7568 

229 

2315 

4240 

kick 100 km 

1592 

690 

2786 

6430 

147 

1928 

3578 

kick 300 km 

393 

172 

666 

1524 

38 

479 

854 

qA=0.01 

385 

287 

783 

1715 

59 

509 

948 

«A=0.1 

1758 

752 

3030 

6862 

183 

2136 

3902 

«A=10 

1810 

834 

3184 

7080 

216 

2200 

4016 

/3=0 

2157 

984 

3707 

8247 

248 

2653 

4649 

/3=0.5 

3067 

677 

4413 

9573 

71 

3402 

5587 

/3=1 

1554 

336 

2249 

4984 

43 

1715 

2843 

Worst 

111 

18 

160 

358 

4 

123 

205 


Table 29. 140 cases with the conservative core-merger criterion, 1000 years, S/N ^ 35 

The same as Table 24 but for 140 cases with the conservative core-merger criterion 


14models 

Inspiral 
S/N^ 35 
1/1000 yrs] 

QNM(0.70) 
S/N3: 35 
1/1000 yrs] 

Quadrature sum of 
Inspiral and 
QNM(0.70) 

S/N> 35 

1/1000 yrs] 

Linear sum of 
Inspiral and 
QNM(0.70) 
S/N^ 35 
1/1000 yrs] 

QNM(0.98) 
S/N;j 35 
1/1000 yrs] 

Quadrature sum of 
Inspiral and 
QNM(0.98) 

S/N^ 35 

1/1000 yrs] 

Linear sum of 
Inspiral and 
QNM(0.98) 
S/N^ 35 
1/1000 yrs] 

our standard 

1373 

4983 

6979 

12584 

1140 

2858 

5385 

IMFdogflat 

1548 

2862 

4879 

9012 

614 

2412 

4321 

IMF:Salpeter 

1030 

714 

1977 

3769 

153 

1263 

2088 

lEFxonst. 

1336 

4820 

6796 

12216 

1095 

2714 

5218 

IEF:e-0® 

1354 

4614 

6574 

11850 

1029 

2694 

5120 

kick 100 km s“l 

1144 

2951 

4591 

9106 

760 

2136 

4157 

kick 300 km s“^ 

313 

1192 

1616 

2749 

284 

643 

1185 

qA=0.01 

170 

244 

480 

1053 

71 

289 

588 

«A=0.1 

1167 

2007 

3695 

7679 

522 

1924 

3901 

«A=10 

1144 

4244 

5915 

10287 

854 

2270 

4246 

/3=0 

1432 

5100 

7136 

12715 

1189 

2921 

5482 

/3=0.5 

1786 

5209 

7792 

14520 

1085 

3206 

6250 

/3=1 

811 

1424 

2596 

5658 

430 

1407 

2853 

Worst 

99 

35 

151 

313 

2 

111 

182 


2030, we can detect gravitational waves from Pop III BH- burst (GRB), the high redshift observation of GRB by Hi-z 
BHs which merged at 2 ^ 9. Thus, we might identify the GUNDAM (Yonetoku et al. 2014) might be possible. Thus, 

peak of Fig. 24. The peak of Fig. 24 depends not only on the merging NS-NS and NS-BH might also be useful for Pop 

binary parameters, but also on the Pop III SFR. Therefore, III binary parameter studying, 
we might get the information of Pop III SFR. 

In this paper, the merging NS-NS and NS-BH are not 
considered because they are negligibly small in number in 
almost all models. However, in some models they are not 
so. Since they are the candidates of the short gamma ray 
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